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Analyses  of  free  vibrations,  forced  vibrations  and  structural  damping 
characteristics  of  sandwich  and  laminated  composites  with  imperfect  interface  and 
delaminations  are  presented.  A refined  high-order  theory  for  the  vibration  of 
sandwich  beams  and  plates  to  specifically  deal  with  a transversely  flexible  core  has 
been  introduced.  The  effect  of  imperfect  viscoelastic  interface  on  the  structural 
damping  of  a sandwich  beam  is  studied.  Analytical  models  of  vibrations  of 
sandwich  and  laminated  plates  with  delaminations  are  developed. 

The  refined  vibration  theory  embodies  a general  systematic  approach  to  the 
analysis  of  sandwich  beam  and  plate  having  high-order  effects  owing  to  the  non- 
linear pattern  of  the  longitudinal  and  transverse  deformation  of  the  thick  core. 
Inertia  effects  due  to  all  the  displacements  are  taken  into  account.  As  such  it 
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improves  on  the  available  elassical  and  shear  deformation  theories.  It  is  used  to 
eorporate  the  imperfect  interface  model  and  vibration  damping  approach  to  exam 
the  interface  layer  effect  on  the  structural  damping  of  the  sandwich  beam.  Results 
show  that  the  imperfect  interface  affects  vibration  damping  significantly. 

Based  on  this  theory  the  free  vibration  behavior  of  a sandwich  plate  with 
a through-width  delamination  is  examined.  The  shear  deformation  effect  and 
compressive  in-plane  force  effect  have  been  taken  into  account  in  the  analysis. 
Two  analytical  models  of  different  deformation  modes  are  discussed.  Fundamental 
frequencies  and  vibration  mode  shapes  are  calculated  and  the  physical  admissible 
model  is  identified.  It  is  found  that  delamination  length,  delamination  location 
and  core  thickness  have  significant  effect  on  vibration  characteristics. 

Finally,  vibration  analysis  of  laminated  composite  beams  with  delaminations 
based  on  a layer-wise  finite  element  model  is  presented.  Parallel  algorithm  design 
for  the  finite  element  analysis  is  discussed  and  efficient  approximation  algorithms 
are  developed.  The  fundamental  natural  frequencies  of  laminated  composites  with 
delamination  are  examined  for  the  effects  of  delamination  length  and  compressive 
in-plane  loading. 
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CHAPTER  1 
INTRODUCTION 


1.1  Background 

Sandwich  construction  has  been  used  in  aeronautical  applications  for  more 
than  40  years,  since  it  has  many  advantages,  for  example,  high  bending  stiffness, 
good  weight  saving,  good  surface  finish,  good  fatigue  properties,  good  thermal  and 
acoustical  insulation,  etc.[l].  Today,  there  is  a renewed  interest  in  using  sandwich 
structures  due  to  the  introduction  of  new  materials  such  as  laminated  composites 
for  the  face  materials  and  various  nonmetallic  honeycomb  and  plastic  foam  core 
materials.  However,  these  new  materials  also  induce  some  new  problems  that  need 
to  be  solved.  Delamination  is  one  of  the  most  frequently  encountered  problems. 
To  use  new  materials  efficiently,  a good  understanding  of  their  structural  dynamic 
behaviors  under  various  loads  is  needed. 

Advanced  fiber-reinforced  composite  materials  such  as  graphite-epoxy 
laminates  offer  several  advantages  for  structures  which  require  a combination  of 
high  strength  and  stiffness  with  low  weight.  With  their  capability  of  being 
tailored  for  a specific  application,  composite  materials  have  been  widely  used  in 
aeronautical  industries  to  replace  metals.  However,  these  materials  are,  in  general, 
more  brittle  than  their  metallic  counterparts.  Also,  in  contrast  to  their  in-plane 
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properties,  transverse  tensile  and  interlaminar  shear  strengths  of  composite 
laminates  are  quite  low.  As  a consequence,  advanced  composite  laminates  are 
susceptible  to  delamination  from  a wide  variety  of  sources. 

Delamination,  also  referred  to  as  interface  cracking  or  debonding  where 
adjacent  laminae  become  unattached  from  one  another,  is  a unique  failure  mode 
in  sandwich  and  laminated  composites.  For  composite  sandwich  panels,  an 
interface  bond  between  the  face  and  the  core  is  relatively  thicker  and  may  be 
weaker  than  that  in  laminated  composites.  In  many  sandwich  structures,  the  bond 
is  not  of  the  perfect  kind  that  can  be  modeled  by  continuity  of  tractions  and 
displacements  across  a discrete  face-core  interface.  Instead,  the  bond  is  effected 
across  a thin  interfacial  zone,  which  has  its  own  mechanical  properties.  Poorly 
bonded  areas  are  probable  locations  for  delaminations  to  occur,  and  these  regions 
of  imperfect  bonding  might  be  attributed  to  foreign  matter  or  nonuniform 
distribution  of  the  adhesive  material.  Clearly,  such  an  interface  may  significantly 
affect  the  integral  dynamic  properties  of  the  sandwich  structures. 

Delamination  may  occur  as  a consequence  of  the  imperfections  in  the 
production  process  such  as  fabrication  stresses,  or  the  effects  of  external  factors 
during  the  operational  life  such  as  environmental  cycling,  handling  damage,  and 
foreign  object  impact  damage.  During  manufacturing,  for  example,  air  pockets 
may  arise  from  layers  not  being  compressed  against  one  another  before  curing 
process.  Once  fabricated,  delaminations  may  develop  from  an  impact  as  minor  as 
a dropped  screwdriver  or  other  tools.  Hence,  the  delamination  can  be  located  in 
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various  places  in  the  sandwich  and  laminate  structures.  Regardless  of  its  location, 
delamination  may  reduce  the  overall  stiffness,  lower  the  load-carrying  capacity  of 
the  composite  panels  and,  thereby,  alter  the  static  and  dynamic  characteristics  of 
the  structural  system  [2,  3].  Sandwich  and  laminate  structures  are  usually  being 
used  in  situations  which  require  them  to  withstand  vibrational  forces.  The 
resonant  frequency  shift  and  degradation  of  dynamic  stiffness  due  to  the 
delamination  may  lead  to  global  structural  failure  at  loads  below  the  design  load. 

The  presence  of  delamination  has  also  been  of  major  concern  in  the  marine 
industry,  where  sandwich  composites  have  attracted  considerable  attention.  For 
example,  highly  advanced  racing  boats  employ  carbon/epoxy  laminates  in  the  hull 
faces  over  low-density  PVC  foam  core  to  reduce  weight  [4].  Their  technology  is 
based  on  wet  lay-up  processes  and  room  temperature  cure  resins  other  than  using 
autoclaves  and  specialized  prepregs  to  achieve  satisfactory  face-core  bonds. 
Although  the  success  has  been  tremendous,  the  marine  industry  is  recognizing  that 
delamination  between  the  core  and  faces  due  to  wave  slamming  on  critical  regions 
of  the  hull  may  present  a problem  [5].  Cyclic  water-generated  loading  or  transient 
wave-slamming  loading  of  high  magnitude  and  short  duration  may  cause  initiation 
and  growth  of  face-core  debonding  that  may  be  disastrous  for  the  load  bearing 
capability  of  the  structure.  The  tendency  for  this  failure  mode  is  expected  to  be 
greatly  enhanced  by  the  presence  of  flaws  or  partial  debonding  that  are  likely  to 
exist  between  the  core  and  face  sheets.  Considerable  vibration  within  the  plate 
components  of  the  structure  can  be  set  up  under  working  conditions  and  there  may 
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have  a severe  effect  upon  the  plate  stability.  It  is,  therefore,  necessary  to 
undertake  vibration  and  dynamic  stability  analyses  and  to  design  the  system 
accordingly. 

Since  dynamic  response  is  governed  by  stiffness,  mass,  and  damping,  the 
effects  of  damage  such  as  delamination  and  crack  on  these  properties  are  of 
interest.  The  use  of  dynamic  testing  to  detect  and  quantify  the  effects  of  internal 
damage  on  modal  parameters  has  been  studied  on  a continuous  but  limited  basis 
over  the  last  twenty  years.  In  recent  years  several  studies  have  been  concerned 
with  damage-induced  changes  in  stiffness  as  found  from  observed  changes  in 
resonant  frequencies  of  composite  specimen  [6,  7].  Other  studies  have  shown  that, 
for  a variety  of  structures  and  loading  conditions,  damping  is  generally  far  more 
sensitive  to  damage  than  stiffness  [8,  9].  More  recently,  the  amount  of  research 
in  this  area  has  increased  due  primarily  to  the  availability  of  inexpensive  multi- 
channel Fourier  analyzers.  Shifts  in  natural  frequencies  have  been  used  to  identify 
damage  generated  by  static  loading  [10,  11],  fatigue  loading  [12,  13],  mechanical 
damage  such  as  saw  cuts  and  local  cracks  [14,  15],  foreign  object  impact  [16-18], 
and  prescribed  delamination  [19-21].  Although  these  work  show  that  system 
dynamic  property  measurements  can  be  used  qualitatively  detect  damage,  there  are 
still  some  major  problems  remaining  and  much  more  effort  will  be  needed  to 
devote  to  this  area.  A better  understanding  of  vibration  and  damping 
characteristics  of  sandwich  and  laminated  structures  with  delamination  will 
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certainly  help  achieve  more  effective  means  of  detecting  the  various  delamination 
damage  using  dynamic  testing. 

1.2  Literature  Survey 

A significant  amount  of  research  has  been  conducted  on  the  vibration  and 
stability  analyses  of  sandwich  and  laminated  beams  and  plates.  Both  analytical 
(closed  form,  Galerkin,  Rayleigh-Ritz  method)  and  numerical  methods  have  been 
used.  A recent  report  by  Leissa  [22]  contains  over  400  references  relative  to  the 
research  work  in  this  area  on  sandwich  and  laminated  plates  and  shells.  A review 
of  the  technical  literature  on  dynamic  behavior  of  laminated  composites  was  given 
by  Bert  [23]  and  on  the  application  of  the  finite-element  method  to  study  the 
vibration  response  of  composite  plates  by  Reddy  [24,  25]. 

1.2.1  Sandwich  Vibration  Theories 

Most  of  the  literature  on  the  vibration  of  sandwich  panels  can  be  separated 
into  two  areas.  The  first  is  to  deal  with  undamped  sandwich  beams,  plates  and 
shells  without  considering  highly  damped  core  effects  on  vibration  response. 
Those  work  are  most  concerned  about  the  calculation  of  vibration  mode  shape  and 
the  evaluation  of  the  critical  buckling  loads  for  various  types  of  instabilities 
including  the  combination  and  interaction  with  face  wrinkling  for  sandwich  panels 
[26].  The  second  aspect  is  to  investigate  the  vibration  response  of  viscoelastically 
damped  sandwich  structures  including  beams,  plates,  rings  and  shells.  Most  of  the 
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research  done  in  this  area  is  to  find  an  effective  analytical  model  of  predicting 
system  loss  factors  and  damped  resonance  frequencies  under  forced  vibration. 

Sandwich  panels  with  traditional  honeycombed  cores  were  studied  by  many 
researchers.  The  basic  approach  used  in  sandwich  static  and  dynamic  problems  is 
based  on  the  assumption  of  an  anti-plane  core,  i.e.,  the  core  with  shear  rigidity 
only,  its  plane  of  section  remains  plan  after  deformation  and  the  height  of  the  core 
is  kept  unchanged  [27,  28].  The  Hexcel  designer  manual  [29]  outlines  design 
procedures  for  such  sandwich  beams.  Ojavo  [30]  assumed  different  deflections  to 
the  two  faces,  but  ignored  the  peeling  stresses.  Ogorkiewicz  [31]  dealt  with  a 
foam  core  by  replacing  it  with  an  ordinary  beam  with  equivalent  properties. 
Weissman-Berman  [32]  developed  an  elastic  foundation  model  to  improve  the 
sandwich  beam  theory  with  simply  supported  boundary  conditions. 

Vibration  problems  of  sandwich  panels  using  the  anti-plane  assumption  were 
solved  either  by  finite  element  method  [33]  or  analytical  approaches  [34,  35].  In 
recent  years,  various  vibration  theories  of  undamped  sandwich  beams  have  been 
developed  by  Ko  [36],  Lu  and  Libove  [37]  and  Bai  and  Sun  [38]  in  the  similar 
basis.  A high-order  theory  used  by  Reddy  [39,  40]  may  be  applied  to  such  panels 
to  overcome  the  deficiencies  encountered  in  the  classical  theories.  These  two 
theories,  however,  assume  that  the  core  is  incompressible  in  the  transverse 
direction,  and  the  boundary  conditions  at  the  same  section  are  identical  for  the  two 
faces.  These  assumption  are  unrealistic  whenever  a sandwich  panel  with  a soft 
compressible  core  is  considered,  particularly  with  concentrated  mechanical 
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impedance  or  partially  distributed  loads.  Those  characteristics  must  be  addressed 
by  a refined  high-order  theory. 

The  damping  effect  has  not  taken  into  account  in  the  above  studies,  and 
only  free  vibration  has  been  considered.  However,  the  damping  capacity  of  the 
composite  sandwiches  is  much  greater  due  to  the  highly  damped  core  materials  and 
their  shear  deformations.  Flexural  vibration  of  viscoelastic  sandwich  beams  has 
long  been  of  interests  since  1960’s.  In  a very  fundamental  paper,  Ross,  Ungar  and 
Kerwin  [41]  treated  the  problem  by  assuming  perfect  interface,  the  same  amount 
of  transverse  displacement  in  each  layer,  damping  only  from  the  shearing  of  the 
viscoelastic  core,  and  simply  supported  boundary  conditions.  DiTaranto  [42] 
extended  this  analysis  to  freely  vibrating  beams  having  any  type  of  boundary 
conditions.  Based  on  DiTaranto ’s  work  sixth-order  differential  equations  of 
motion  were  derived  in  terms  of  transverse  displacement  by  Mead  and  Markus 

[43] ,  and  this  theory  was  correlated  with  experimental  data  by  Lu  and  Douglas 

[44] .  Yan  and  Dowell  [45]  used  the  principle  of  virtual  work  to  obtain  a set  of 
five  differential  governing  equations.  A more  accurate  theory  including  the  inertia 
effects  of  transverse,  longitudinal  and  rotary  motions  was  carried  out  by  Rao  and 
Nakra  [46]. 

Although  there  are  some  published  papers  dealing  with  the  delamination 
problems  of  sandwich  panels,  to  the  author’s  best  knowledge,  none  of  them 
consider  the  dynamic  response.  The  buckling  and  postbuckling  of  a delaminated 
sandwich  beam  with  laminated  faces  and  honeycomb  core  were  studied  by  Hwu 
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and  Hu  [47].  Debonding  failure  of  sandwieh  eomposite  materials  has  been 
investigated  analytieally  and  experimentally  by  Carlsson,  Sendlein  and  Merry  [48]. 
They  employed  laminated  beam  shear  deformation  theory  to  calculate  effective 
bending  properties  of  delamination  region  and  fracture  toughness.  The  study 
presented  here  is  trying  to  conduct  some  basic  vibration  analyses  of  delaminated 
sandwich  beams  and  plates  with  various  core  and  face  materials. 

1.2.2  Interface  Effects  of  Sandwich  Panels 

Existing  theories  of  vibration  damping  of  sandwich  structures  have 
invariably  assumed  perfect  bonding  between  faces  and  core  as  well  as  constant 
transverse  displacement  through  the  thickness  of  the  beams.  Subsequently  the 
plane  of  sections  remains  plane  and  the  transverse  normal  (peeling)  stresses  are 
neglected.  On  the  other  hand,  the  assumption  of  negligible  normal  stress  of  the 
core  (or  constant  transverse  displacement  through  the  core)  is  not  applicable  when 
the  core  is  made  of  low  strength  materials  such  as  foams  and  aramid  honeycomb. 
In  reality,  however,  for  flexible  core  materials  a more  complex  state  exists 
between  faces  and  core.  The  assumptions  of  equal  amount  of  transverse 
displacement  and  perfect  bonding  may  not  be  suitable  to  describe  the  problem 
accurately. 

The  existence  of  an  interface  layer  or  adhesive  layer  between  the  face  and 
the  core  on  the  physical  behavior  of  the  sandwich  composite  materials  has  been 
recognized  during  recent  years.  The  interface  has  been  described  as  either  a layer 
of  specified  thickness  and  of  elastic  constants  different  from  those  of  the  face  and 
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the  core,  or  a very  thin  interfacial  zone  of  unspecified  thickness  where  tractions 
are  continuous  but  the  displacements  may  not.  Advanced  racing  boats,  for 
example,  employ  carbon/epoxy  laminates  in  the  hull  skin  over  low  density  PVC 
foam  where  interfacial  layers  are  considered  to  be  effective  in  energy  absorption 
[49].  Fazio,  Hussein  and  Ha  [50]  introduced  a nonrigid  bonding  effect  into  the 
elasticity  solutions  for  sandwich  beam-columns  under  transverse  static  loading. 
Imperfect  interface  conditions  have  been  quantitatively  defined  as  linear  relations 
between  interface  displacement  discontinuities  and  associated  traction  components. 
Hussein  [51]  presented  analytical  solutions  of  sandwich  plates  assembled  with 
nonrigid  adhesive  layers.  Their  results  have  shown  that  the  bonding  stiffness,  up 
to  a certain  level,  has  a strong  effect  on  the  deflections  and  stresses  of  beam- 
columns  and  plates.  Such  adhesive  layer  effects  on  buckling  load  was  also 
extended  to  consider  the  case  of  viscoelastic  interface  between  viscoelastic  cross- 
ply  composite  faces  and  the  elastic  core  by  Kim  and  Hong  [52].  In  the  present 
study,  we  will  examine  the  effect  of  viscoelastic  interface  on  the  structural 
damping  of  sandwich  beams. 

1.2.3  Laminated  Structure  Vibration  Theories 

Kapania  and  Raciti  [53,  54]  have  given  a more  recent  survey  on  vibration 
response  of  symmetric  and  unsymmetric  laminated  beams  and  plates.  We  hereby 
give  some  representative  research  work  in  this  field.  Reddy  and  Kuppusamy  [55] 
performed  a vibration  analysis  of  laminated  plates  using  a three-dimensional 
solution.  Bhimaraddi  and  Stephens  [56]  presented  results  of  vibrations  of  thick 
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laminated  plates  using  a high-order  theory.  The  effect  of  transverse  shear  on  the 
cylindrical  bending,  vibration,  and  buckling  of  laminated  plates  was  presented  by 
Stein  and  Jegley  [57].  Effects  of  delamination  on  vibration  and  stability 
characteristics  of  laminated  composites  have  not  been  considered  until  recent 
years. 

The  effect  of  delamination  on  buckling  load,  postbuckling  deformation  and 
delamination  growth  with  various  geometrical  parameters,  loading  conditions, 
material  properties  and  boundary  conditions  have  been  studied  extensively.  Most 
studies  on  the  subject  of  delamination  buckling  and  growth  were  based  on  the 
static  postbuckling  solutions  of  homogeneous  or  laminated  plates  [58,  59].  The 
energy  release  rate  is  used  in  the  analysis  and  the  Griffith-type  fracture  criterion 
with  a constant  specific  fracture  energy  is  assumed  in  these  studies.  These  include 
the  differentiation  of  the  total  strain  energy  with  respect  to  the  delamination  length 
by  Chai  et  al.  [60]  and  Chow  and  Ngan  [61],  a finite  element  method  to  compute 
the  crack  closure  integral  by  Whitcomb  [62],  and  the  path  independent  J integral 
by  Yin  and  Wang  [63].  Dynamic  analysis  of  the  growth  for  a thin,  homogeneous, 
strip  delamination  in  a thick,  axially  loaded  base  plate  was  studied  recently  [64- 
66]. 

A more  general  theory  for  modeling  of  delamination  in  composite  laminates, 
based  on  a layer-wise  displacement  field,  was  presented  by  Barbero  and  Readdy 
[67].  They  suggested  that  the  same  displacement  distributions  in  individual  layer 
is  capable  of  representing  displacement  discontinuity  conditions  at  interfaces 
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between  layers.  It  demonstrated  that  the  layer-wise  laminate  theory  provides  an 
adequate  framework  for  the  delamination  analysis  of  laminated  composites.  It 
must  be  noted,  however,  that  the  computational  cost  of  the  analysis  makes  it 
unattractive  for  prediction  of  delamination  buckling  behavior  in  engineering 
practice  when  the  number  of  layer  is  large.  Based  on  a 1-D  layer-wise  theory,  a 
finite  element  model  was  developed  by  Lee  et  al.  [68]  to  examine  the  buckling  of 
axially  loaded  composite  beams  with  multiple  through-width  delaminations.  A 
parallel  study  by  Suemasu  [69,  70]  shows  common  numerical  trends  that  are 
similar  to  those  observed  by  Lee  et  al.  Kutlu  and  Chang  [71]  also  studied  the 
compression  failure  of  laminated  composites  containing  multiple  through-width 
delaminations.  The  effect  of  elliptical  delamination  of  composite  plates  on  elastic 
and  buckling  behavior  was  investigated  by  Peck  and  Springer  [72].  However,  only 
a few  investigations  of  delamination  effects  on  vibration  characteristics  have  been 
conducted. 

Kulkarni  and  Frederick  [73]  were  among  the  first  to  study  the  vibration 
characteristics  of  laminated  structures  with  delaminations.  They  considered  the 
particular  case  of  a circular  cylindrical  shell  with  a circumferentially  symmetric 
delamination  crack  of  small  length  at  the  middle  surface.  In  their  study,  the  effect 
of  the  delamination  was  taken  into  account  as  a reduced  bending  rigidity  of  the 
shell  over  the  delamination  zone,  obtained  by  replacing  the  moment  of  inertia  of 
the  delaminated  layers.  However,  as  show  in  some  later  studies,  such  a modeling 
grossly  overestimates  the  reduction  in  the  bending  stiffness  caused  by  the 
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delamination.  Using  a similar  model,  Ramkumar  et  al.  [74]  formulated  the 
vibration  of  a delaminated  beam  based  on  the  Timoshenko  beam  theory.  The 
coupling  effect  of  the  oscillatory  axial  forces  in  the  upper  and  lower  delaminated 
layers  was  not  taken  into  account  in  the  analysis  which  results  in  consistently 
much  lower  natural  frequencies  as  compared  to  the  experimental  data.  This 
consistency  exists  even  for  cases  with  short  delamination.  Wang  et  al.  [75]  then 
studied  the  vibration  of  isotropic  beams  with  a through-width  delamination, 
analyzing  the  flawed  beam  as  four  different  Euler  beams  that  are  joined  together, 
with  appropriate  boundary  and  matching  conditions  imposed  to  obtain  the  response 
of  the  beam  as  whole.  The  coupling  effect  of  axial  forces  in  split  layers  was 
considered  in  the  analysis.  Their  analysis  examined  various  delamination  length 
and  locations,  and  their  effects  on  the  beam  natural  frequencies  and  mode  shapes. 
It  was  found  that  for  the  delaminated  beam  with  a short  and  close-to-midplane 
delamination,  the  calculated  fundamental  frequency  is  close  to  the  experimental 
results.  Their  analysis,  however,  leads  to  modes  of  deformation  of  the  delaminated 
beam  which  are  not  physically  admissible  in  the  case  of  off-midplane 
delaminations  because  the  sublaminates  above  and  below  the  plane  of  delamination 
were  implicitly  assumed  to  be  free  to  deflect  and  thus  have  different  transverse 
displacements  in  the  formulation. 

For  vibration  relative  to  the  natural  (unloaded)  referential  state  in  which 
there  is  no  opening  between  the  upper  and  lower  delaminated  layers,  geometric 
incompatibility  (i.e.  overlap  between  upper  and  lower  sublaminates)  will  occur  if 
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the  mode  shapes  of  the  upper  and  lower  sublaminates  are  different.  A constraint 
of  identical  transverse  displacements  of  upper  and  lower  sublaminates  have  been 
imposed  by  Mujumdar  and  Suryanarayan  [76]  in  their  proposed  analytical  model. 
Their  analysis  is  based  on  a one-dimensional  Euler  beam  analysis  in  such  a 
"constrained  mode"  can  take  into  account  delaminations  located  at  an  arbitrary 
position  in  both  the  spanwise  and  thicknesswise  directions.  Numerical  results  for 
frequencies  were  obtained  by  using  this  "constrained  mode"  for  various  boundary 
conditions  and  mode  numbers  for  a certain  range  of  values  of  the  delamination 
size,  spanwise  location  and  thicknesswise  location  of  the  delamination. 

Vibration  response  of  delaminated  composite  beam-plates  relative  to  the 
postbuckling  state  has  also  been  studied.  Yin  and  Jane  [77]  considered  the 
vibration  behavior  of  a compressed  clamped-clamped  beam-plate  with  a 
symmetrically  located  delamination  in  the  buckled  states.  The  natural  frequencies 
of  the  delaminated  plate  at  postbuckled  states  were  compared  to  the  natural 
frequencies  of  a perfect  plate  for  various  delamination  sizes  and  depths.  The 
effect  of  delamination  in  a cantilevered  composite  beam  on  natural  frequencies  and 
mode  shapes  is  evaluated  both  analytically  and  experimentally  by  Shen  and  Grady 
[78]  recently.  For  vibrations  relative  to  the  postbuckling  state,  the  constraint  of 
identical  transverse  displacements  of  upper  and  lower  sublaminates  is  not 
necessary  because  there  is  an  opening  between  the  upper  and  lower  sublaminates. 
A free  mode  formulation  as  used  in  Wang’s  work  [75]  is  appropriate.  Chen  and 
Goggin  [79]  used  shear  deformation  theory  to  determine  the  transverse  shear  effect 
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on  natural  frequencies  and  mode  shapes  in  compressed  delaminated  plates  at  the 
postbuckled  state. 

It  can  be  seen  that  most  of  the  research  on  the  vibration  response  of 
laminated  composites  with  delaminations  is  restricted  to  the  symmetric  delaminated 
cases.  Symmetric  delaminated  laminates  are  simpler  to  analyze  since  they  do  not 
exhibit  bending-extension  coupling  that  is  characteristic  of  unsymmetric  laminates. 
When  symmetric  laminates  are  delaminated  during  manufacturing  or  operation, 
unsymmetric  laminates  are  usually  formed.  As  a result,  the  coupling  effect  may 
be  important  on  the  vibration  characteristics  of  delaminated  plate.  Besides  there 
are  some  other  reasons  for  designing  structural  elements  with  such  coupling 
effects.  As  suggested  by  Minich  and  Chamis  [80],  for  example,  it  would  be 
possible  to  select  coupling  effects  that  will  provide  built-in  self-damping 
mechanisms  when  subjected  to  dynamic  responses.  Also  the  coupling  effects  are 
used  in  aeroelastic  tailoring  to  control  mode  shapes  and  frequencies  of  high-aspect- 
ratio  lifting  surfaces  [81,  82].  In  the  present  study,  the  bending-extension 
coupling  effects  on  the  vibration  response  of  delaminated  laminates  will  be 
included. 

1.3  Objective  and  Outline  of  Present  Study 

Accurate  knowledge  of  dynamic  behavior  of  sandwich  and  laminated 
composites  is  essential  for  reliable  and  lightweight  structural  design,  particularly 
in  the  high-speed  transportation  industries.  The  presence  of  imperfect  interface 
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and  delamination  not  only  affect  the  strength  and  integrity  of  the  structure  by 
contributing  to  the  final  failure,  but  also  cause  a reduction  in  the  stiffness,  thus 
affecting  its  vibration  and  stability  characteristics.  As  was  discussed  in  the 
previous  section,  there  have  been  very  few  research  work  in  the  subject  of 
vibration  characteristics  of  laminated  structures  with  delaminations,  and  no  study 
toward  the  vibration  of  sandwich  panels  with  delaminations  has  been  conducted 
to  date.  In  the  research  area  of  sandwich  constructions,  the  viscoelastic  face-core 
interface  effect  has  also  not  been  examined  in  the  dynamic  response  region. 

The  main  objective  of  this  research  is  to  develop  analytical  models  for  the 
analysis  of  dynamic  responses  of  sandwich  and  laminated  composites  with 
imperfect  interface  and  delamination.  Emphasis  is  placed  on  developing  refined 
beam  and  plate  vibration  theories  which  is  specially  suitable  for  the  sandwich 
construction  with  compressible  and  flexible  core  materials.  Another  problem  of 
interest  is  to  examine  viscoelasic  interface  effects  on  structural  damping  of 
sandwich  panels  and  free  vibration  behavior  of  delaminated  sandwich  and 
laminated  composites.  A finite  element  model  based  on  layer-wise  theory  is  also 
implemented  in  conjunction  with  parallel  algorithm  design  for  finite  element 
analysis. 

In  the  present  study,  refined  beam  and  plate  theories  of  sandwich  vibration 
damping  are  derived  via  variation  of  the  Hamilton  principle.  The  nonlinear 
displacement  field  of  the  core  in  both  longitudinal  and  vertical  directions  through 
the  thickness  are  assumed  based  on  the  elasticity  solution.  The  face  layers  are 
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considered  as  the  ordinary  beams,  acting  longitudinally  only  and  interconnected 
with  the  core  through  equilibrium  and  compatibility  at  the  interface  layers.  The 
boundary  conditions  are  given  at  each  face  layers  and  core  respectively, 
accomplished  by  either  prescribing  the  stress  resultants  or  the  corresponding  time 
derivatives  of  the  displacement  components.  Imperfect  interface  between  face  and 
core  of  a sandwich  panel  is  quantitatively  defined  as  the  longitudinal 
displacements  discontinuities  at  the  interface  layers  being  proportional  to  the 
associated  shear  stresses.  Vibration  analysis  of  delaminated  sandwich  panels  is 
performed  using  this  refined  theory. 

The  organization  of  the  dissertation  is  as  follows:  in  the  next  chapter  we 
derive  a refined  vibration  theory  of  sandwich  beam  using  Hamilton  principle  and 
perform  the  verification  and  parametric  studies.  An  interface  model  is  introduced 
and  viscoelastic  face-core  interface  effects  on  vibration  damping  of  sandwich 
beams  are  examined  in  Chapter  3.  Chapter  4 extends  the  refined  vibration  beam 
theory  into  a plate  theory.  In  Chapter  5 vibration  of  delaminated  sandwich  plates 
is  studied  based  on  the  refined  theory.  Vibration  analysis  of  laminate  beams  is 
carried  out  in  Chapter  6 using  a finite  element  model  based  on  the  layer-wise 
theory  with  a parallel  computation  consideration.  Finally,  conclusions  drawn  from 
this  study  and  recommendations  for  future  research  are  presented  in  Chapter  7. 


CHAPTER  2 

A REFINED  VIBRATION  THEORY  OF  SANDWICH  BEAMS 


The  effects  of  more  flexible  core  materials  on  the  dynamic  response  and 
damping  properties  of  sandwich  beams  are  studied  by  introducing  a refined 
vibration  theory  of  sandwich  beams.  The  theory  of  the  core  is  based  on  the 
elasticity  solution  of  more  accurate  displacement  assumptions,  including  nonlinear 
transverse  deformation.  The  face  layers  are  considered  as  ordinary  beams  with 
axial  and  bending  resistance. 

With  the  viscoelastic  properties  of  the  core  material  which  is  complex  and 
is  function  of  frequency  for  a given  temperature,  the  linear  equations  of  motion 
describing  the  vibration  of  the  sandwich  finite  beam  are  derived  based  on  the 
Hamilton  principle.  In  addition  to  transverse  inertia  effects,  longitudinal 

translatory  and  rotary  inertia  effects  have  been  included.  Consistent  boundary 
conditions  are  accomplished  by  either  prescribing  the  stress  resultants  or  the 
corresponding  time  derivatives  of  the  displacement  components.  Typical  cases  of 
a simply  supported  beam  under  the  harmonic  excitation  are  numerically 
investigated.  The  numerical  results  of  mechanical  impedance  examples  are  well 
agreed  with  experimental  data  in  the  literature.  The  effects  of  core  thickness  and 
stiffness  on  resonance  frequencies  are  studied  parametrically. 
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2.1  Basic  Assumptions 

Consider  a three  layer  sandwich  beam  of  length  / and  unit  width,  subjected 
to  a transverse  dynamic  load  q{x,t).  It  has  face  layers  of  thickness,  /i,  and  /ij,  and 
a core  of  thickness  h shown  in  Figure  2.1.  The  face  layers  are  the  elastic  material 
with  high  Young’s  moduli  of  E,  and  Ej.  The  flexible  core  is  linearly  viscoelastic 
with  a Young’s  modulus  and  a shear  modulus  of  where  E^  and  are 
expressed  by  complex  numbers  under  the  harmonic  excitation. 


Z 


Figure  2.1  Sandwich  Beam  Construction 
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The  following  assumptions  underlie  the  present  development: 

• The  elastic  face  layers  are  considered  as  ordinary  beams  with  axial  and 
bending  resistance. 

• The  viscoelastic  core  carries  null  longitudinal  stresses,  but  takes  nonlinear 
displacement  fields  in  both  x and  z directions. 

• The  perfect  interface,  i.e.,  displacement  continuity  is  hold  at  the  face-core 
interface. 

Based  on  the  assumptions,  the  following  displacement  fields  are  given  at  the 


top  and  bottom  face  layers: 


u'(x,z,t)  = u,(x,t)  - z. 

dx 


dwpcj) 


(2.1) 


w‘(x,z,t)  = w^(x,t) 


dwp,t) 


(2.2) 


w%x,z,t)  = w^{x,t) 
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For  the  core: 


u‘^{x,z,t)  = u^(x,t)  - z + zei|;(x,0 


5% 


dx 


+ 


(2.3) 


1 dx  6 


2 dx 


where  superscripts  t,  b and  c denote  top,  bottom  face  layers  and  the  core,  Ui(x,t) 
and  Wj(x,t)  (i  = 0,  1,  2)  are  the  longitudinal  and  vertical  displacements  of  the 
centroid  of  the  core  and  face  layers,  e = E^G^,  is  a constant  and  Zj  = Z - h - hj/2, 
Z2  = Z + h2/2  are  measured  from  the  centroid  of  the  faces  respectively,  as  shown 
in  Figure  2.2.  Deformation  pattern  and  sign  convention  for  stresses,  resultants  and 
displacements  also  appear  in  Figure  2.2. 

The  core  displacement  field  is  derived  from  the  three-dimensional  elasticity 
solution  with  only  an  assumption  of  null  in-plane  stress  (assumption  2).  The 
assumption  used  here  follows  those  encountered  in  the  classical  small  elastic 
deflection  theory  of  sandwich  panels  where  the  core  has  shear  rigidity  only  and  its 
resistance  to  in-plane  longitudinal  stresses  is  negligible.  However,  the  transverse 
normal  stresses  (peeling  stresses)  are  considered  in  the  core  as  well  as  in  the 
interface  layers.  The  cross-section  plane  of  the  core  may  distort  and  its  height 
may  change.  It  can  be  seen  from  Eqn  (2.3)  that  the  transverse  displacement  of  the 


core  can  take  parabolic  variations  with  the  core  thickness. 
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(a)  Coordinate  System 


(b)  Sign  Convention 

I I I I I I I I I I I I I ■ I I L 


(c)  Deformation  Pattern 

Figure  2.2  Coordinate  and  Deformation  of  the  Sandwich  Beam 
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2.2  Governing  Equations  of  Motion 


The  kinematic  relations  for  the  face  layers,  using  small  deformation  theory, 
can  be  written  as 


t 


8 


X 


(2.4) 


b 


8 


JC 


(2.5) 


The  kinematic  relations  for  the  core  are: 


e!  = (|)(x,/)  - z 


y1-  = ^ 


d ij;  (x,  t) 
dx 


(2.6) 


The  compatibility  conditions  of  deformations  between  the  faces  and  the  core 
are  described  by  continuity  of  displacements  across  a face-core  interface. 


At  z=h  (upper  interface),  we  have 


(2.7) 

(2.8) 


and  at  z = 0 (lower  interface). 
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(2.9) 

(2.10) 


Substituting  Eqns.  (2.1)-(2.3)  into  above  equations,  Uj{x,t),  U2(x,t) 

and  W2(x,t)  can  be  expressed  in  terms  of  Wj(x,t),  W2(x,t),  UQ(x,t)  and  \|;(a:,0, 


(X,  t)  = Wg  (X,  t) 


(2.11) 


2 dx 


0 


w + /z(j) 


dij; 
2 dx 


(2.12) 


(2.13) 


(2.14) 


A set  of  force  and  moment  resultants  are  defined  for  top  and  bottom  faces 


separately: 
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2 


(2.15) 


2 


where  i = t,  b correspond  to  y = 1,2  respectively. 

Then  the  constitutive  relations  for  the  elastic  face  layers  can  be  written  as 
follows. 


where  and  £,/,  are  the 
For  the  viscoelastic 


, du 

OX 


K - 


dx 


m!  = -EJ 


11 


d^w 

dx^ 


dw. 


h ^ -r  ^ 2 


K - 


dx 


axial  and  the  flexural  rigidities 


core,  the  constitutive  equations 


(2.16) 


of  the  face  layers, 
are  in  the  form. 


T 


c 


xz 


(2.17) 


where  and  are  complex  constants. 
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The  governing  equations  of  motion  and  the  boundary  conditions  are  derived 
via  variational  calculus,  the  Hamilton  principle,  which  requires: 


where  T is  the  kinetic  energy;  U is  the  total  potential  energy;  t,,  define  the  time 
interval. 

The  total  kinetic  energy  of  the  sandwich  beam  can  be  expressed  as: 


where  p,,  pj  and  p are  the  density  of  the  upper  face,  lower  face  and  the  core, 
respectively,  and  ( ) denotes  a first  derivative  with  respect  to  time. 

The  total  potential  energy,  U,  consists  of  the  contributions  from  both  face 
layers  Uj  and  U2,  core  shear  strain  U^,  core  normal  U4,  and  the  external  potential 
term  V due  to  applied  loading. 


(2.18) 


(2.19) 


u - a ^ u,  ^ u,  ^ u,  ^ V 

12  3 4 


(2.20) 
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(2.21) 


F = / ' qw^dx  (2.22) 

Jo 

where  q = q{x,t)  is  transverse  distributed  dynamic  loading  at  surface  of  top  face 
layer. 

It  is  noted  that  T and  U can  be  further  expressed  in  terms  of  four  basic 
unknown  functions — w^(x,t),  W2{x,t),  Uq{xJ)  and  i|;(A:,r).  Taking  variation  in  the 
Hamilton’s  equation  with  respect  to  w,,  m>2,  Uq  and  ij;,  four  equations  of  motion  and 
consistent  boundary  conditions  are  obtained  after  integration  by  parts  in  x and  t, 
and  some  algebraic  manipulations  and  substitutions.  It  is  noted  that  homogeneous 
conditions  for  the  displacement  and  the  velocity  at  ? = are  assumed.  The  four 
governing  equations  can  be  written  in  the  matrix  form: 


[dM\{u] . [di:]{c/}  . {/=•) 


(2.23) 
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where 


{U}  = {w„  Wj,  Uo, 


{F}  = {q,  0,  0,  0}^ 


[dM\  — the  mass  matrix  differential  operators 
[3^  — the  stiffness  matrix  differential  operators 
It  is  noted  that  [dK],  {U}  and  {U}  are  complex  valued.  In  appendix  A the  mass 
matrix  differential  operators  [dM]  and  stiffness  matrix  differential  operators  {dK\ 
are  given  in  terms  of  the  four  independent  variables,  w,,  W2,  Mq  t);.  It  can  be 
seen  that  inertia  effects  due  to  all  the  assumed  displacements  are  taken  into 
account. 

The  boundary  conditions,  at  the  end  of  the  beam,  are  given, 
at  the  top  face  layer  ( x = 0 and  x = I ): 


= 0 or  u = 0 


(2.24) 


dw 


= Q or 


0 


(2.25) 


dx 


+ 


dx  12  dx  213  4 


h dw^  3hdw^ 


(2.26) 


10  5a:  20  dx  45  dx^ 


or  w,  = 0 
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at  the  bottom  face  layer  ( x = 0 and  x = I ): 


=0  or  = 0 


(2.27) 


= 0 or 


dx 


0 


(2.28) 


0JC 


'2^  ^2 

- p 

• • 

Seh 

12  dx 

3 

24 

3h^x 

2h  ^^2 

13A'  a'i{r 

40  dx 

15  dx 

720  ax' 

\ 


= 0 


/ 


(2.29) 


or  ^ 


and  at  the  core  { x — Q and  x — I ): 


ij;  = 0 


or 


h 5\j; 


12  dx 
/ 


ph 


3E 


\2E 


C L 


Wj  + + 


/ 


\ 


\ 


7e 


/ 


aij; 

10  dx 


\ 


1 _h_d^  \3h^  d^i|j 

4 ax  ’ 15  ax'  120  ax'  S40  ax' 


\ 


= 0 


/ 


(2.30) 


It  is  noted  that  the  boundary  conditions  of  the  core  at  each  end  must  be  satisfied 
at  any  point  along  the  height  of  the  core.  It  can  also  be  observed  that  boundary 
conditions  that  differ  from  one  face  layer  to  the  other  at  the  same  section  are 
allowed  here.  In  Appendix  A all  boundary  conditions  expressed  in  terms  of  four 


independent  functions  are  given. 
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2.3  Solution  of  Forced  Harmonic  Motion 


Consider  a sandwich  beam  under  the  harmonically  varying  loading, 


q{x,t)  = Q(x)e"^‘ 


(2.31) 


where  w is  the  forcing  frequency.  In  general,  as  the  frequency  of  q(x,t)  varies, 
both  the  amplitude  and  mode  of  vibration  vary.  We  assumed  here  that  the 
response  due  to  the  applied  force  is  harmonic  and  vibrating  at  the  excitation 
frequency. 

Let  {U(x,t)j  be 


w,(x) 


{f/(x,0} 


(2.32) 


Thus  the  partial  differential  operators  in  Eqn  (2.23)  change  to  ordinary 
differential  operators,  and  the  governing  equations,  Eqn  (2.23)  becomes: 


(-(o^[M]  . = {/}  (2.33) 

where  {/(x)}  = (Q(x),  0,  0,  0}'^. 

The  set  of  simultaneous  ordinary  differential  equations  described  by  Eqn 
(2.34)  can  be  solved  using  known  procedures  from  mathematics  [85],  in 
conjunction  with  the  boundary  conditions  as  indicated  in  Eqns  (2.24)  to  (2.31). 
Since  the  motion  is  now  harmonic,  it  is  legitimate  to  admit  viscoelastic  damping 
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into  the  core  by  putting  = G'(l  + ina)  + *tie)-  The  manner  in 

which  the  loss  factor,  tiq  and  iig,  and  the  real  part  of  the  shear  modulus,  G'  and 
E'  vary  with  temperature  and  frequency  is  quite  different  for  various  viscoelastic 
materials. 

It  should  be  pointed  out  that  due  to  the  introduction  of  viscoelastic  damping 
material  in  the  system,  all  physical  quantities  of  the  system  are,  in  general, 
complex.  Although  the  purely  real  solutions  for  the  displacements,  for  example, 
are  not  a general  possibility,  there  are  nonetheless  exceptions  [43]. 

Since  one  of  the  primary  objectives  of  employing  the  damping  composite 
structures  is  to  attenuate  the  vibratory  responses  at  the  resonant  conditions,  to  find 
the  resonance  frequencies  of  such  a system  is  of  great  importance.  However,  the 
usual  procedure  used  in  vibration  analysis  to  find  the  frequencies  can  hardly  be 
followed  for  such  a composite  system  because  the  characteristic  equation  for  each 
mode  is  difficult  to  obtain.  In  addition,  due  to  the  linear  superposition  application 
and  the  nature  of  the  frequency-dependent  material  properties,  to  obtain  the 
resonance  frequencies  for  such  a damped  system  is  not  easy.  It  is  believed  that 
using  the  forced  vibration  analysis  is  probably  the  most  straightforward  procedure 
to  determine  the  resonance  frequencies.  The  frequencies  can  be  located  when  the 
slopes  of  responses  between  two  frequencies  change  signs  in  the  neighborhood  of 
resonant  conditions.  Nonetheless,  a method  of  determining  resonance  frequencies 
of  damped  beams  will  be  presented  in  the  following  when  the  refined  vibration 
theory  is  used. 
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Consider  a three-layered  damped  sandwich  beam  with  simply  supported 
boundary  conditions.  In  this  case,  all  boundary  conditions  are  satisfied  and  the 
series  form  solution  may  be  expressed  as: 
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where  constants  w,„,  W2„,  Uq„,  i|;„  (n  = 1,  2,...)  are  the  amplitude  coefficients  and 
n is  the  wave  number. 

Substituting  Eqn  (2-34)  into  Eqn  (2-23),  after  some  manipulation,  and 
utilizing  the  orthogonality  relations  of  the  sine  and  cosine  functions,  four  linear 
frequency  dependent  algebraic  equations  governing  the  displacement  components 
tVi„,  Wjn,  Uq„,  and  i|;„  can  be  found.  Then,  the  transverse  displacement  component 
of  the  upper  face  may  yield. 


iV(fa)) 

Z)((o) 


(2.35) 


The  resonance  frequencies  can  be  obtained  by  solving  the  characteristic  equation 


Z)(co)  = 0 


(2.36) 
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The  solution  of  this  equation  yields  four  frequeneies  for  a given  mode.  Each  of 
the  eigenvalues  in  general  is  expected  to  be  complex  because  the  coefficients  of 
Eqn  (2.36)  are  complex.  Three  family  modes  of  vibration  are  predominantly 
flexural,  extensional  and  thickness  shear  types.  For  each  mode  number,  the 
resonance  frequency  and  modal  damping  (modal  loss  factor)  can  be  calculated. 

For  the  damped  system,  Eqn  (2.35)  can  be  reduced  to  such  a form  that  the 
real  and  imaginary  quantities  are  present  only  in  the  denominator.  Following  the 
usual  procedure  in  vibration  analysis  for  small  damping,  the  resonance  frequencies 
for  each  mode  can  be  determined  by  solving  the  polynomial  of  the  denominator 
with  the  imaginary  part  neglected.  However,  for  the  case  in  which  there  is  a 
considerable  amount  of  damping  in  the  system,  the  aforementioned  procedure  to 
locate  the  resonance  frequencies  may  not  be  applied  because  the  resonance 
frequencies  may  not  be  close  to  the  solutions  obtained  based  on  the  small  damping 
assumption  of  vibration  analysis.  For  generality,  a procedure  is  proposed  to  find 
the  frequencies  based  on  the  definition  of  the  resonant  condition  for  damped 
systems. 

The  resonant  condition  for  the  damped  system  is  defined  as  the  frequency 
for  which  the  response  is  maximum.  For  each  mode,  the  resonance  frequency  may 
be  obtained  by  differentiating  the  amplitude  of  Eqn  (2.35)  with  respect  to  the 
frequency  parameter  and  setting  it  equal  to  zero.  Thus 
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or 
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a 0)'  = 0 


(2.37) 


This  frequency  equation  turns  out  to  be  a 15th-order  polynomial  with  real 
coefficients.  Because  of  lengthy  expression,  the  explicit  form  of  this  polynomial 
is  not  given  here.  Solving  this  polynomial,  the  resonance  frequency  for  each  mode 
can  be  determined. 


2.4  Verification  Problems 

To  verify  the  new  refined  vibration  theory  of  damped  sandwich  beams,  a 
free-free  damped  sandwich  beam  subjected  to  a concentrated  time  harmonic  load 
with  circular  frequency  of  excitation  o at  the  mid-span  was  examined  to  correlate 
both  experimental  data  and  analytical  solutions  published  by  Lu  and  Douglas  [44]. 
In  their  experiment,  two  damped  three-layer  laminates  were  constructed  to 
compare  with  the  analytical  solutions  given  by  Mead  and  Markus  [43]. 

Specimen  1 incorporated  a thin  0.004-inch  (0.01  cm)  acrylic  base 
viscoelastic  material,  with  mass  density  1.03  x lO"'*  lb-s^/in'‘  (1.124  x 10'^  g-sVcm''), 
sandwiched  between  two  identical  steel  machined  flats,  1 /4-inch  thick  (0.635  cm), 
1-inch  wide  (2.54  cm),  and  24  3/16-inches  long  (61.456  cm).  The  machined  flats 
helped  ensure  a uniform  cross  section  for  the  thin  damping  core.  The  viscoelastic 
material  used  is  a 3M  Company  product,  two-sided  adhesive  tape  designated  as 
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SJ2015X.  The  shear  storage  modulus  and  shear  loss  modulus  of  the  viscoelastic 
layer  were  obtained  for  the  frequency  region  of  50  Hz  to  1000  Hz  from 
measurements  made  by  the  manufacturer  on  a Melabs  rheometer,  an  instrument 
that  sinusoidally  excited  small  rectangular  specimens  in  a low  amplitude  shear 
mode.  The  complex  shear  modulus  at  124°  (51°)  could  be  functionally 
approximated  by  the  expressions 


InG'  = 0.5791n—  . 1.136  (G'inpsi) 

2n 


InG"  = 0.601  In—  . 1.144  (G"inpsi) 

2ti: 


(2.38) 


where  the  shear  modulus  G*  = G(1  + iti)  = G'  + iG"  (1  psi  = 6895  Pascal). 
Examination  of  the  temperature  and  frequency  dependence  of  the  complex  shear 
modulus  showed  that  this  material  approximated  a class  A thermorheologically 
simple  material.  By  using  this  assumption,  the  complex  shear  modulus  was 
extended  to  3000  Hz. 

Specimen  2 incorporated  a 1/4-inch  (0.635  cm)  thick  neoprene  rubber 
damping  core  with  mass  density  9.35  x 10'^  Ib-sVin'*  (1.02  x 10'^  g-sVcm'*) 
sandwiched  between  two  steel  faces,  each  1/4-inch  (0.635  cm)  thick,  2 inches 
(5.108  cm)  wide,  and  18  1/8-inches  (46.038  cm)  long.  Comparison  of  the  result 
of  specimen  1 with  the  thicker  damping  core  in  this  specimen  permitted  an 
evaluation  of  the  analytical  model  over  a wide  range  of  elastic  layer  to  viscoelastic 
layer  thickness  ratios.  With  values  given  in  reference  [86],  the  complex  shear 
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modulus  at  77°  (25°)  was  found  to  functionally  approximate  the  expressions, 


(0 


InG'  = 0.026 In—  . 4.754  (G'  in  psi) 

2tc 


InG"  = O.nOln—  - 2.705  (G"  in  psi) 

2;t 


(2.39) 


for  the  frequency  range  of  50  Hz  to  10000  Hz  (1  psi  = 6895  Pascal). 

The  beams  were  excited  by  an  electromechanical  vibration  generator 
wrapped  around  a mechanical  impedance  head  which  was  attached  to  the  beam  at 
mid-span  by  an  area-reducing  cone  to  concentrate  the  dynamic  load.  The  beams 
were  suspended  in  air  by  lightweight  elastic  shock  cords  which  resulted  in  a 
system  natural  frequency  well  below  5 Hz.  The  force  and  acceleration  signals 
were  analyzed  by  2 Hz  bandwidth  tracking  filters,  and  the  20  Hz  to  3000  Hz 
spectrum  was  swept  slowly  with  a sweep  oseillator  to  obtain  a eontinuous 
speetrum.  The  acceleration  signal  was  integrated  to  obtain  velocity.  The 
oscillator  sweep  rate  was  slow  enough  to  allow  the  strueture  and  eleetronics  to 
aecurately  record  impedance.  The  beams  were  tested  in  a temperature  controlled 
chamber  after  thermal  equilibrium  was  reached. 

The  analytical  solutions  of  those  mechanical  impedance  examples  are 
obtained  for  the  free-free  damped  sandwich  beams  by  using  the  present  refined 
vibration  theory.  The  terminology  "mechanical  impedance"  used  here  is  the  ratio 
of  force  causing  motion  to  the  velocity  at  a given  point  in  a vibrating  system.  The 
coneentrated  load  applied  at  the  mid-span  is  expressed  in  terms  of  infinite  series. 
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The  eight  free-free  boundary  conditions,  first  equation  of  Eqns  (2.24)  through 
(2.31)  are  used,  and  the  complete  solution  can  be  obtained  as: 
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(2.40) 


where  the  a/s  are  the  eight  roots  of  the  characteristic  equation  with  complex 
coefficients,  the  Afs  are  the  eight  unknown  complex  constants  to  be  determined 
from  boundary  conditions.  The  eight  unknown  complex  constants  can  be 
computed  by  any  standard  matrix  inversion  technique,  such  as  the  Gauss- Jordan 
elimination  method.  Once  the  transverse  displacements,  w^(x,t),  is  found,  the 
driving  point  mechanical  impedances  can  be  determined.  Since  the  core  used  here 
is  relatively  thin,  the  results  show  that  w^ix.t)  and  W2{x,t)  are  almost  the  same. 

Figures  2.3  and  2.4  show  the  comparison  for  the  mechanical  impedance 
versus  frequency  between  experimental  data  and  present  analytical  results.  These 
figures  are  plots  of  absolute  values  of  the  driving  point  mechanical  impedances 
versus  frequency  in  Hertz.  The  two  specimens  selected  are  in  very  different 
conditions.  The  first  has  a very  thin  relatively  highly  damped  viscoelastic  core 
(SJ2015X)  measured  at  124°F  (51°),  while  the  second  has  a relatively  thick  low 
loss  tangent  viscoelastic  core  (neoprene  rubber)  measured  at  77°  (25°).  In  both 
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Figure  2.3  Driving  Point  Mechanical  Impedance  of  a Damped  Sandwich  Beam  (Specimen  1) 
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Figure  2.4  Driving  Point  Mechanical  Impedance  of  a Damped  Sandwich  Beam  (Specimen  2) 
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cases,  not  only  are  the  damped  resonance  frequencies  well  predicted  by  the  refined 
theory  but  also  the  resonant  amplitudes.  It  is  noted  that  although  Mead  and 
Markus  theory  [43]  is  well  agreed  with  the  experiment  data,  the  present  theory 
improves  the  high  mode  results. 

2.5  Parametric  Study 

For  a given  damped  system,  it  is  anticipated  that  the  vibrational  behavior 
is  complicated  by  the  involvement  of  many  parameters.  It  depends  on  many 
geometrical  and  material  parameters  of  the  structural  system  including  the 
modulus,  width,  length,  and  thickness  of  the  upper  and  lower  face  layer  structures; 
the  thickness,  shear  modulus  and  loss  factor  of  the  damped  core  material;  the  half 
wavelength  of  the  mode  which  in  turn  depends  on  the  mode  number  and  the 
boundary  conditions.  Moreover,  the  damping  material  properties  are  functions  of 
frequency  as  well  as  temperature  which  has  made  the  design  of  engineering 
structures  with  damping  for  a specific  application  even  more  difficult.  Therefore, 
for  an  optimized  design  of  any  particular  damped  structural  system,  a complete 
parametric  study  is  required. 

The  resonance  frequency  behavior  of  a simply  supported  sandwich  beam  is 
studied  parametrically.  The  results  are  presented  for  a symmetric  construction 
which  consist  of  two  identical  faces  and  a soft  core.  The  non-dimensional  n-th 
mode  frequency  parameter  used  in  the  figures  is  defined  as 
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2 2 9hJ 

Q„  = w (1  + ir\  ) 

” ” " E,/, 


(2.41) 


The  fundamental  mode  resonates  at  the  frequency  w,  and  has  a loss  factor  of  ti,. 
Unless  mentioned  otherwise,  the  loss  factors  of  the  core  are  assumed  as  tiq  = 0.2 


and 


0.05. 


The  effect  of  the  core  thickness  on  the  resonant  frequency  is  presented  in 
Figure  2.5.  It  can  be  seen  that  the  resonant  frequency  decreases  with  increasing 
the  core  thickness,  and  Mead  and  Markus  theory  overestimates  the  resonant 
frequency  especially  for  the  thicker  core. 

Figure  2.6  shows  the  effect  of  the  modulus  of  elasticity  of  the  core  in  the 
vertical  direction  on  the  fundamental  resonant  frequency.  In  the  Mead  and  Markus 
theory,  = 0,  meaning  that  = <»,  so  that  it  overestimates  the  resonance 
frequency.  Such  discrepancy  may  be  greater  if  the  core  material  is  getting  softer. 

In  Figure  2.7,  the  effect  of  the  modulus  of  elasticity  of  the  core  in  the 
vertical  direction  on  the  resonance  frequency  is  examined  in  the  other  way  in 
which  the  frequency  is  normalized  with  respect  to  the  ordinary  beams.  The 
frequency  of  a Euler-Bernoulli  beam  with  the  same  mass  and  flexural  rigidity  is 
denoted  by  Oq,  and  that  of  a Timoshenko  beam  which  has  shear  rigidity  and 
rotational  inertia  is  denoted  by  Wj.  The  damping  effects  are  neglected  here  and 
first  and  second  mode  are  considered.  It  can  be  seen  that  shear  rigidity  and 
rotational  inertia  have  less  effects  on  the  second  mode  than  the  first  mode.  It  is 
noted  that  an  increase  in  the  core  modulus  E^  is  always  associated  the  with  an 
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Figure  2.5  Variation  of  First  Resonance  Frequency  with  Core  Thickness 
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Figure  2.6  Effect  of  Core  Modulus  on  the  Resonance  Frequency 
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Figure  2.7  Shear  and  Rotational  Effects  on  Resonance  Frequencies 
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increase  the  core  density  for  the  foam  or  honeycomb  core  materials.  Thus,  the 
natural  frequency  will  increase  with  the  core  modulus.  The  frequency  ratio 
relative  to  the  homogeneous  isotropic  beam  may  offer  a help  for  the  design  of  low 
density  sandwich  construction  in  the  dynamic  environment. 

According  to  the  definition  of  Mead  and  Markus  [43],  the  shear  parameter 
is  defined  as 
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which  is  an  important  parameter  to  represent  the  core  shear  effect.  Figure  2.8 
examines  the  variation  of  fundamental  resonance  frequency  with  the  shear 
parameter  s in  different  values  of  face/core  thickness  ratio.  It  can  be  seen  that  the 
resonance  frequency,  and  hence  the  stiffness  of  the  beam,  changes  more  rapidly 
with  s when  the  face/core  thickness  ratio  is  smaller.  This  may  be  significant  to 
the  performance  of  a sandwich  beam  under  changing  temperature  because  of 
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Figure  2.8  Variation  of  Resonance  Frequencies  with  Shear  Parameter 


CHAPTER  3 

EFFECTS  OF  VISCOELASTIC  INTERFACE  OF  SANDWICH  BEAMS 


The  effects  of  viscoelastic  interface  layer  (or  adhesive  layer)  on  the 
dynamic  response  and  structural  damping  of  sandwich  structures  are  studied  by 
employing  a refined  vibration  theory  of  sandwich  beams  as  discussed  in  Chapter 
2.  Imperfect  interface  conditions  between  faces  and  core  are  defined  as  linear 
relations  between  interface  displacements  discontinuities  and  associated  shear 
stresses.  With  the  viscoelastic  properties  of  the  interface  and  the  core  material 
which  are  complex  and  is  function  of  frequency  for  a given  temperature,  the  linear 
equations  of  motion  describing  the  vibration  of  the  sandwich  finite  beam  are 
derived  based  on  the  Hamilton  principle. 

Three  methods  to  evaluate  system  loss  factors  (structural  damping)  are 
discussed.  In  those  methods  the  overall  damping  of  the  structure  is  calculated 
from  the  independent  properties  of  the  constituent  materials  that  form  the  sandwich 
structure.  The  effects  of  storage  modulus  and  loss  factor  of  interface  layers  on  the 
structural  damping  and  resonance  frequencies  are  numerically  investigated  for 
simply  supported  beams  under  the  harmonic  excitation.  Results  demonstrate  that 
the  influence  of  viscoelastic  interface  layer  on  damping  is  significant  and  there 
may  exist  optimal  adhesive  bonding  stiffness  to  achieve  the  maximum  sandwich 
structural  damping. 
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3.1  Viscoelastic  Interface  Model 

Consider  a three  layer  sandwich  beam  of  length  I and  unit  width  with  two 
interface  layers  between  upper,  lower  face  layers  and  core  as  shown  in  Figure  3.1. 
It  has  face  layers  of  thickness,  hj  and  h2,  and  a core  of  thickness  h.  The  face 
layers  are  the  elastic  material  with  high  Young’s  moduli  of  E,  and  Ej.  The 
flexible  core  is  linearly  viscoelastic  with  a Young’s  modulus  E^  and  a shear 
modulus  of  Gj.  where  E^  and  are  expressed  by  complex  numbers  under  the 
harmonic  excitation. 

The  classical  compatibility  conditions  of  deformations  between  the  face  and 
core  assumed  the  displacements  and  associated  stress  are  continuous  across  the 
interface  layer.  This  will  be  referred  to  as  perfect  interface  conditions  and  every 
other  kind  of  interface  conditions  will  be  called  imperfect.  Since  most  of  energy 
dissipation  in  the  viscoelastic  damping  mechanism  is  contributed  by  the  shear 
deformation,  the  effect  of  transverse  normal  deformation  at  the  interface  layer  on 
structural  damping  is  neglected,  and  the  continuity  conditions  of  transverse 
displacement  at  the  interface  layer  remain.  The  imperfect  interface  conditions  are 
described  by  the  longitudinal  displacement  discontinuity  as  described  in  the 
following  assumption 

• The  longitudinal  displacement  discontinuity  (slip)  is  proportional  to  the 


shear  stress  at  the  face-core  interface  which  is  viscoelastic. 
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Figure  3.1  A Sandwich  Beam  with  Interface  Layers 
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To  study  the  structural  damping  problem  due  to  the  viscoelastic  interface 
layers  and  core,  the  complex  modulus  approach  is  used  and  the  compatibility 
conditions  can  be  written  as  the  following 
at  z = h (upper  interface) 

= k'  {u‘  - u ‘^)  (3.1) 


w"  = w' 


at  z = 0 (lower  interface) 

T2  = k'  - w *) 


(3.2) 


(3.3) 


w = w * (3.4) 

where  x and  x 2 are  the  shear  stress  at  the  upper  and  lower  interface  layers 
respectively,  and  k = k(w)(l  + is  the  frequency-dependent  complex  shear 

modulus  of  interface  layers.  Small  values  of  k((n)  represent  the  weak  bonding 

between  the  faces  and  the  core.  When  k — ► °o,  it  approaches  perfect  bonding 

interface  that  is  usually  assumed  in  the  sandwich  theory. 

3.2  Governing  Equations  of  Motion 

A set  of  force  and  moment  resultants  are  defined  for  top  and  bottom  faces 


separately: 
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(3.5) 


where  i = t,  b correspond  to  y = 1,2  respectively. 

The  constitutive  relations  for  the  elastic  face  layers  and  for  the  viscoelastic 
core  are  exactly  the  same  form  as  we  discussed  in  Chapter  2,  equations  (2.16)  and 
(2.17)  where  frequency-dependent  complex  shear  and  Young’s  modulus  of  the 
core,  G(.,  E^,  are  used. 

Eliminating  Uo(x,t),  vVo(x,t),  and  4>(jc,t)  from  Eqns  (3.1)  through  (3.4),  an 
independent  compatibility  equation  is  obtained: 
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= 0 


The  governing  equations  of  motion  and  the  boundary  conditions  are  derived 
via  variational  calculus,  the  Hamilton  principle,  which  requires: 
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[T  - U] 


dt  = 0 


(3.7) 


where  T is  the  kinetic  energy;  U is  the  total  potential  energy;  define  the  time 


interval. 
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The  total  kinetic  energy  of  the  sandwich  beam  can  be  expressed  as: 


dw 

OX 


+ 


(3.8) 


where  p,,  pj  and  p are  the  density  of  the  upper  face,  lower  face  and  the  core, 
respectively,  and  the  ( ) denotes  a first  derivative  with  respect  to  time. 

The  total  potential  energy,  U,  consists  of  the  contributions  from  both  face 
layers  C/,  and  U2,  core  shear  strain  U^,  core  normal  U^,  and  potential  terms  V,  due 
to  applied  loading.  The  effects  of  imperfect  bonding,  (p  is  considered  by 
introducing  a Lagrange  multiplier  X with  Eqn  (3.6)  into  the  total  potential  energy 
contributions.  Thus, 


U = ^ ^ ^ ^ (P  ^ V 


(3.9) 
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(3.10) 


cp  = l\‘ 
Jo 


h,  + h dw,  h.■^  h dw. 

k\u,  -u,  . — ^ -) 

' ^ 2 ax  2 dx 


{k  he^2E  )yU  .k—^ 

' 12  ax' 


dx 


V = I qw  dx 
Jo 


(3.11) 


where  q = q{x,t)  is  transverse  distributed  dynamic  loading  at  surface  of  top  face 
layer. 

It  is  noted  that  T and  U can  be  further  expressed  in  terms  of  five  basic 
unknown  functions — w,(x,t),  W2{x,t),  u^(x,t),  U2(x,t)  and  v|;(x,t).  Taking  variation 
in  Hamilton’s  equation  with  respect  to  w^,  W2,  u^,  U2,  lit  and  X,  six  equations  of 
motion  and  boundary  conditions  are  obtained  after  integration  by  parts  in  x and  t, 
and  some  algebraic  manipulations  and  substitutions.  It  is  noted  that  homogeneous 
conditions  for  the  displacement  and  the  velocity  at  t = t,,  ?2  assumed.  The 
physical  meaning  of  Lagrange  multiplier  X is  the  face-core  interface  shear  force. 
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The  number  of  governing  equations  is  reduced  from  six  to  five  by  eliminating  the 
A.  The  governing  equations  can  be  written  in  the  matrix  form: 

[aM]{?/} . [a^]{t/}  = {F}  (3-12) 


where 

{U}  = {w^,  W2,  u^,  U2, 

{F}  = {q,  0,  0,  0,  0}^ 

{dM\  — the  mass  matrix  differential  operator 

[0^  — the  stiffness  matrix  differential  operator 

• • 

It  is  noted  that  [dK\,  {U}  and  {U}  are  complex  valued. 

The  boundary  conditions,  at  the  end  of  the  beam,  are  given, 
at  the  top  face  layer  ( x = 0 and  x = I ): 

(3.13) 

(3.14) 
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(3.15) 


54 


at  the  bottom  face  layer  ( x = 0 and  x = I ): 
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(3.16) 

(3.17) 


(3.18) 


and  at  the  core  ( x = 0 and  x = I ): 
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It  should  be  noticed  that  the  boundary  conditions  of  the  core  at  each  end  must  be 
satisfied  at  any  point  along  the  height  of  the  core.  Boundary  conditions  that  differ 
from  one  face  layer  to  the  other  at  the  same  section  are  allowed  here.  It  is  also 
noted  that  all  boundary  conditions,  traction  type  or  displacement  type,  can  also  be 
expressed  in  terms  of  five  basic  functions  (see  Appendix  B). 

The  eight  simply  supported  boundary  conditions  at  each  end  of  the  beam 
will  be  a combination  of  these  boundary  conditions  as  the  following 
at  the  top  face  layer,  (x  = 0 and  x = 1) 

Nj  - 0 (3.20) 

Mj  - 0 (3.21) 

w,  = 0 (3.22) 


at  the  bottom  face  layer, 

(3.23) 

(3.24) 

(3.25) 


K - 0 
m/  = 0 

w,  = 0 


at  the  core. 


i|/  = 0 (<_.  = 0) 


(3.26) 
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Similarly  the  fixed  end  and  free  end  boundary  conditions  can  be  obtained 
in  a combination  of  those  boundary  conditions  given  in  Eqns  (3.13)  through 
(3.19). 

3.3  Calculation  of  Loss  Factors 


Structural  damping  refers  to  a structure’s  or  structural  component’s  capacity 
for  dissipating  energy  or  to  its  capacity  for  removing  from  a structural  vibration 
some  of  the  energy  associated  with  that  vibration.  This  removed  energy  may  be 
converted  directly  into  heat,  transferred  to  connected  structures  or  ambient  media. 
The  mechanism  of  structural  damping  is  a complicated  physical  phenomenon.  In 
addition  to  controlling  the  amplitude  of  resonant  vibrations,  damping  also  has 
other  effects  such  as  increasing  wave  attenuation,  decreasing  sound  transmission 
through  structures,  reducing  structural  fatigue,  etc. 

Usually,  a convenient  measure  of  the  damping  of  any  dynamical  system  is 
its  non-dimensional  loss  factor  which  can  be  defined  by  the  energy  dissipated  per 
cycle  of  harmonic  motion  in  a particular  mode  divided  by  the  maximum  potential 
energy  stored  in  the  course  of  the  cycle  multiplied  by  2n, 


U 

n = 

2%W 


(3.27) 


where  U is  the  energy  dissipated  per  cycle  or  the  energy  that  has  to  be  supplied 
to  sustain  the  steady  state  vibration,  and  W is  the  maximum  energy  associated  with 
vibration.  However,  there  are  many  quantitative  measures  of  structural  damping 
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ability  based  on  the  effect  it  has  on  the  response  of  the  system.  Most  of  the 
measures  of  damping  that  are  used  currently,  such  as  resonance  bluntness,  half- 
power bandwidth,  logarithmic  decrement,  percentage  of  critical  damping  etc,  are 
defined  based  on  linear  single  degree  of  freedom  systems  with  frequency 
independent  parameters.  Of  the  many  different  measures  of  damping  used  in 
resonant  analysis,  loss  factor  and  quantities  closely  related  to  it  are  still  the  most 
generally  applicable.  On  the  other  hand,  all  the  different  measures  of  damping, 
under  the  assumption  of  low  damping,  are  used  interchangeably. 

A lot  of  research  work  on  vibration  damping  only  deal  with  structural 
damping,  also  called  hysteretic  damping,  which  was  originally  postulated  as  a basis 
for  describing  the  internal  damping  properties  of  solids.  It  is  showed  that  the 
structural  damping  can  be  considered  adequately  by  using  a complex  modulus,  E' 
+ iE  , where  the  real  part  is  called  the  storage  modulus  and  the  imaginary  part  the 
loss  modulus.  Loss  factor,  a measure  of  damping,  is  the  ratio  of  the  imaginary 
part  to  the  real  part.  In  the  case  of  structural  damping,  both  E'  and  E"  are 
constants.  Although  the  complex  modulus  is  only  defined  for  harmonic  excitation, 
the  complex  modulus  approach  can  also  be  used  for  the  analysis  of  a system 
subjected  to  arbitrary  excitation  because  any  real  signal  can  be  represented  as  a 
Fourier  sum. 

However,  structural  damping  is  not  a physically  realizable  energy 
dissipation  mechanism  in  a continuous  composite  system  where  it  can  also  be 
modeled  by  using  complex  modulus  approach  based  on  the  correspondence 
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principle  of  viscoelasticity  [85].  It  is  known  that  complex  modulus  which  does  not 
vary  with  frequency  leads  to  noncausal  system  response  (response  antieipates 
input)  [86],  and  a finite  energy  dissipation  rate  at  zero  frequency.  In  the  ease  of 
viscoelastic  damping,  both  E'  and  E"  of  viscoelastic  materials  are  funetions  of 
frequency,  instead  of  eonstants. 

Viscoelastic  damping  in  a continuous  composite  system  ean  also  be  analyzed 
by  using  the  complex  modulus  approach.  According  to  the  elastie-viscoelastic 
correspondenee  principle,  the  elastic  constants  of  the  materials  are  simply  replaced 
by  the  corresponding  frequency-dependent  viseoelastic  counterparts.  In  a 
eonsequence,  the  system  loss  factor  can  be  ealculated  from  the  resulting  eomplex 
valued  system  of  equations  in  the  following  three  approaches. 

3.3.1  Complex  Eigenvalue  Approaeh 

Consider  a sandwich  beam  under  the  harmonieally  varying  loading, 

q{x,t)  - Q{x)e"^‘  (3-28) 

where  w is  the  forcing  frequency.  In  general,  as  the  frequeney  of  q{x,t)  varies, 
both  the  amplitude  and  mode  of  vibration  vary.  We  assumed  here  that  the 
response  due  to  the  applied  force  is  harmonic  and  vibrating  at  the  excitation 


frequency. 
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Let  {U(x,t)}  be 


{C/(x,0}  = 


w,(x) 

W2(X) 

u (x) 
U^ix) 


>e'“' 


(3.29) 


Thus  the  partial  differential  operators  change  to  ordinary  differential 
operators,  and  the  governing  equations,  Eqn  (3.12)  becomes: 

(- = {/■}  (3.30) 


where  {f(x)}  = {0,  Q(x),  0,  0,  0}'^. 

The  eigenvalue  problem  is  then 

[/:][C/J  . 0)/[An[C/J  (3.31) 


where  g)„  and  [[/„]  are  the  n-th  complex  eigenvalue  and  eigenvector. 

For  a undamped  beam,  the  Raleigh’s  principle  (Raleigh’s  quotient)  can  be 
written  as 


^ [ujmu] 

[Uf[M\[U„] 


(3.32) 


and  for  a damped  structural  system,  the  complex  Raleigh’s  principle  can  be 


extended  to 
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{[K^.i[KmUJ  - (0^(1  Wri„)[M][t/J  (3.33) 

where  [^']  is  the  sum  of  the  stiffness  matrix  of  the  elastic  elements,  [^J,  and  the 
real  part  of  the  stiffness  of  the  viscoelastic  elements,  [^rvI-  [^"]  is  the  imaginary 
part  of  stiffness  matrix  of  the  viscoelastic  elements  which  is  equal  to  The 

modal  loss  factor  of  the  n-th  mode  of  the  structure  is  denoted  by  r\„.  Then,  Eqn 
(3.33)  yields 


o>l[M\[U„]  - [K'][UJ 


(3.34) 


where  t|  is  the  material  loss  factor  of  the  core.  Thus, 

^2  [Un][KnU„] 

[U^]  [M]  [C/J 

[U:][K'W„]  [U:][KJ[U„] 

[U„"][M][U„]  [U„^][M\[U„] 

Then 

^[U„'][KJ[U]  F(v) 

^ = n-^ 


(3.35) 


(3.36) 


where  the  strain  energy  of  the  core  associated  with  a given  mode  is  V„(v),  and  the 
strain  energy  of  the  entire  elastic  system  associated  with  the  same  mode  is  V„. 
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It  should  be  noted  here  that  the  loss  factors  of  the  structural  system  only  describe 
a part  of  the  effectiveness  of  the  damped  structure. 

In  solving  the  vibration  problem,  the  solutions  can  be  determined  in  terms 
of  damped  normal  modes.  Both  the  eigenvalues  and  the  associated  eigenvectors 
will  of  course  be  complex  in  nature.  Therefore,  the  disadvantages  of  this  approach 
stem  from  the  inherent  nature  of  the  involvement  of  the  damping  material. 
Namely,  the  computational  effort  is  quite  expensive,  and  the  properties  of  real 
viscoelastic  materials  are  frequency  dependent  instead  of  being  constant  for  storage 
moduli. 

3.3.2  Direct  Frequency  Response  Approach 

In  the  direct  frequency  response  method,  a forced  vibration  solving 
technique  is  used.  The  overall  damping  of  the  system  is  calculated  from  a 
theoretical  response  spectrum  which  is  generated  over  a range  of  frequencies  being 
considered.  Since  the  complex  modulus  approach  is  used,  the  Eqn  (3.30)  is 
complex.  For  each  frequency  of  excitation,  the  displacement  per  unit  applied  force 
is  calculated  by  solving  the  system  of  complex-valued  simultaneous  equations,  and 
the  response  function  is  thereby  generated.  Loss  factor  is  calculated  from  the  real 
part  of  generated  response  spectrum  as  shown  in  Figure  3.2,  or  by  using  the  half- 
power-bandwidth technique.  These  two  methods  are  used  in  general  purpose  finite 
elemenet  analysis  package  for  the  layered  damped  structures,  for  example,  the 
NASTRAN  package.  The  advantages  are  that  these  programs  can  analyze 
structures  with  any  complexity  and  that  they  are  well  developed  and  managed. 
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The  disadvantage  is  that  the  computation  cost  is  expensive.  But  the  major 
difference  between  these  damped  systems  and  other  elastic  undamped  laminated 
structures  analyzed  by  the  finite  element  direct  frequency  response  method  is  that 
the  results  of  these  damped  problems  can  be  obtained  for  only  one  frequency  at 
a time  because  the  shear  modulus  and  the  loss  factors  of  the  viscoelastic  materials 
are  functions  of  frequency  for  a given  temperature.  As  a result,  NASTRAN’s 
standard  normal  mode  analysis  can  not  be  used  to  determine  the  resonance 
frequencies  and  mode  shapes  of  such  systems.  The  analytical  solution  is  thus 
necessary  and  more  accurate  analysis  will  be  helpful  in  understanding  the 
structural  damping  mechanism. 

The  system  loss  factor  can  be  obtained  more  directly  by  using  direct 
frequency  response  method.  As  we  discussed  in  section  2.4  of  Chapter  2,  the 
resonance  frequencies  can  be  obtained  by  solving  the  characteristic  equation 

Z)((o)  = 0 (3.37) 

The  solution  of  normal  damped  mode  can  also  be  calculated.  Then  loss  factor  can 
be  calculated  by  using  the  energy  ratio  method.  After  the  resonant  frequency  and 
the  corresponding  displaced  shape  at  resonance  have  been  calculated,  the  loss 
factor  of  the  beam  system  can  be  calculated  from 


Displacement  (Real  Part) 
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Figure  3.2  Real  Part  of  the  Response  Spectrum 


{c/'fMlt/) 


where  {U*}  is  the  complex  conjugate  of  {U}  and 


(3.38) 


m = [K']  . i[K"] 


(3.39) 
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3.3.3  Modal  Strain  Energy  Approach 

The  concept  of  modal  strain  energy  approach  was  formulated  from  the 
damping  energy  and  the  elastic  strain  energy  of  the  damped  structures.  The 
approach  is  essentially  an  extension  of  Raleigh’s  principle  into  the  complex 
domain.  An  approximate  solution  for  the  loss  factor,  Eqn  (3.36),  can  be  calculated 
by  approximating  the  complex  eigenvector  [f/J  by  the  real  vector  [[/„'],  which  is 
calculated  from  purely  elastic  analysis,  i.e.,  by  suppressing  the  imaginary  part  of 
[K"].  It  is  under  this  major  assumption  that  the  modal  loss  factors  of  the  damped 
structural  system  can  be  easily  determined.  Thus,  this  technique  is  valid  only  for 
systems  with  relatively  small  levels  of  damping  where  the  mode  shapes  and 
frequencies  of  the  damped  and  the  undamped  structure  are  similar.  Assuming  that 
modal  loss  factors  are  a function  of  the  strain  energy  distributions  in  undamped 
eigenmodes  (which  implies  relatively  moderate  material  damping  and  ignores 
intermodal  coupling  effects),  a modal  loss  factor  can  be  defined  for  a given 
eigenmode  of  a damped  system. 

It  was  found  that  loss  factor  estimates  based  on  weighted  ratios  of  strain 
energies  in  the  damped  and  undamped  materials,  with  strain  energy  terms  defined 
as  functions  of  forced  structural  deflection  shape,  agreed  with  loss  factors  based 
on  bandwidths  of  discrete  resonant  response  peaks.  This  discovery  made  loss 
factor  estimation  possible  for  a wide  range  of  modes  exhibited  by  finite  element 
idealizations  of  practical  structures,  where  the  stiffness  matrix  and  the 
corresponding  nodal  displacements  are  real  since  only  the  real  eigenvalue  problem 
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is  solved.  As  a result,  extensive  applieations  of  the  modal  strain  energy  numerical 
method  can  be  easily  used  for  damped  structure  analysis  and  particularly  for 
concept  evaluations  and  for  design  applications.  This  makes  damping  design 
analysis  more  accessible  and  economical  since  undamped  eigenmode  computation 
is  much  less  difficult  conceptually  and  much  less  expensive  than  the  forced 
damped  frequency  response  computation. 


3.4  Numerical  Examples 


Vibration  damping  of  simply  supported  sandwich  beams  is  numerically 
studied.  In  this  case,  all  boundary  conditions  are  satisfied  and  the  solution  is 
analytical  and  may  be  expressed  as: 
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where  w,„,  W2„,  u,n,  Ujn,  v|in  are  the  amplitude  coefficients  and  n is  the  wave 
number. 


To  examine  the  effects  of  viscoelastic  interfacial  bonding,  parametric 


studies  are  carried  out  for  the  variations  of  non-dimensional  resonant  frequency 
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and  loss  factor  of  the  beam  system.  The  effect  of  the  storage  shear  modulus  of  the 
interface  layer  for  a simply  supported  symmetrical  beam  on  the  damping  of 
sandwich  structural  system,  ris,  is  presented  in  Figure  3.3.  It  can  be  seen  that  the 
interface  storage  shear  modulus  (adhesive  bonding  stiffness)  has  a strong  effect, 
up  to  a certain  level,  on  the  beam  loss  factor.  Beyond  this  level,  the  usual 
assumption  of  perfect  bonding  used  in  classical  theories  is  quite  acceptable.  This 
level  to  constitute  "perfect"  interface  is  given  in  terms  of  the  ratio  of  the  adhesive 
stiffness  to  the  core  stiffness.  There  may  exist  the  difference  of  the  level  between 
the  static  response  and  dynamic  response.  On  the  other  hand,  there  may  exist  an 
optimal  value  of  interfacial  shear  modulus  to  achieve  the  maximum  structural 
damping. 

Figure  3.4  shows  the  effect  of  the  loss  factor  of  interface  layer  on  the 
structural  damping.  It  is  shown  that  the  damping  of  the  interface  generally 
enhances  the  structural  damping. 

In  Figure  3.5,  the  effect  of  the  core  thickness  on  the  system  loss  factor  is 
examined.  It  demonstrates  that  increasing  core  thickness  will  generally  enhance 
the  sandwich  system  damping.  The  variation  is  also  dependent  upon  the  interface 
shear  modulus. 

In  Figure  3.6,  the  non-dimensional  square  resonance  frequency  of  the 
sandwich  beam  is  plotted  against  the  aspect  ratio  of  the  core.  It  can  be  seen  that 
the  resonance  frequency  decreases  with  increasing  the  core  thickness,  and  the 
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Figure  3.3  Effect  of  Storage  Shear  Modulus  of  Interface  On  Sandwich 

Beam  Structural  Loss  Factors 


1.0 


h//  = 0.15 


0.001 


h//  = 0.10 


0 


0.2 


0.4 


0.6 


0.8 


1.0 


ri,  / ri 

'k  'c 


Figure  3.4  Effect  of  Loss  Factor  of  Interface  on  Sandwich  Beam 

Structural  Loss  Factors 


69 


0 5 10  15  20  25 


h/hi 


Figure  3.5  Effect  of  Core  Thickness  on  Sandwich  Beam  Structural  Loss  Factors 
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Figure  3.6  Effect  of  Interface  Shear  Modulus  on  First  Resonance  Frequency 
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storage  shear  modulus  of  the  interface  layer  also  has  some  effects  on  the  resonance 
frequencies. 

It  should  be  mentioned  that  as  damping  (loss  factor)  of  the  structure  is 
increased,  the  stiffness  of  the  structure,  in  general,  will  be  decreased.  Therefore, 
by  only  improving  damping  will  not  reduce  the  dynamic  responses  such  as 
displacements,  stresses  and  noise  level  of  the  structure.  It  is  a trade-off  situation 
between  damping  and  stiffness.  The  sandwich  dynamic  stiffness  is  observed  to 
increase  with  increasing  the  storage  modulus  of  the  interface  in  Figure  3.7.  There 
is  also  a certain  value  of  the  interface  modulus  to  keep  the  sandwich  stiffness  close 
to  the  "perfect  bonding"  value.  Generally,  a proper  value  of  interface  shear 
modulus  for  the  trade-off  between  enhancement  in  damping  and  decrease  in 


stiffness  can  be  obtained  in  most  situations. 
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Figure  3.7  Effect  of  Storage  Shear  Modulus  of  Interface  on  Sandwich 

Beam  Structural  Stiffness 


CHAPTER  4 

A REFINED  VIBRATION  THEORY  OF  SANDWICH  PLATES 


In  this  chapter,  the  refined  vibration  theory  of  sandwich  beams  is  extended 
to  the  sandwich  plate  theory.  By  means  of  variational  methods,  the  basic 
equations  of  vibratory  bending  of  sandwich  plates  with  laminated  faces  and 
orthotropic  core  are  developed.  Taken  into  consideration  are  the  effects  of 
flexural  and  membrane  energies  in  the  faces,  transverse  shear  and  normal  strain 
energies  in  the  core,  and  rotatory,  translator  and  transverse  inertias  in  both  core 
and  faces.  These  effects  are  found  to  be  of  importance  only  at  very  high 
frequencies  for  homogeneous  beams  and  plates.  These  effects  are,  however,  of 
considerable  importance  at  the  lower  frequencies  for  the  case  of  sandwich  beams 
and  plates. 


4.1  Basic  Formulations 

Consider  a sandwich  plate  with  anisotropic  composite  laminate  faces  and 
an  orthotropic  honeycomb  or  an  isotropic  foam  core.  The  geometry  of  the 
sandwich  plate  considered  is  shown  in  Figure  4.1.  The  thickness  of  the  core  is 
much  larger  than  that  of  the  face,  and  the  transverse  modulus  of  the  core  normal 
to  the  plate  is  much  lower  than  that  of  the  face.  The  flexible  core  is  linearly 
viscoelastic  with  a Young’s  modulus  E^  in  the  z direction  and  a shear  modulus  of 
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and  where  E^,  G^^  and  G^^  are  expressed  in  terms  of  complex  numbers  under 
the  harmonic  excitation. 


Z 


Figure  4.1  Geometry  of  a Sandwich  Plate 

In  order  to  determine  the  response  of  a sandwich  plate  due  to  harmonic 
excitation,  the  assumptions  are  identical  to  those  for  sandwich  beams,  as  given  in 
sections  2.1  and  3.1.  The  classical  laminate  plate  theory  is  used  to  deal  with  the 
face  sheets  as  the  displacement  assumption  gives  in  the  following. 
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For  the  top  face  layer: 


,,  , , , dw(x,y,i) 

u {x,y,z,t)  = u^{x,y,t)  - z, 

ox 


dw^(x,y,t) 

v‘{x,y,z,t)  = v^{x,y,t)  - z, 

dy 


w\x,y,z,t)  = w,(x,y,0 


For  the  bottom  face  layer: 


4,  , , , dw^{x,y,t) 

u {x,y,z,t)  = u^{x,y,t)  - z^ 

OX 


b,  ^ .X  dwp,y,t) 

V (x,y,z,t)  = v^{x,y,t)  - Zj ; 

dy 

w\x,y,z,t)  = w^{x,y,t) 


(4.1) 


(4.2) 


where  Zj  = Z - h - hj/2,  Z2  = Z + h2/2  are  measured  from  the  centroid  of  the  faces 
respectively,  as  shown  in  Figure  4.1.  The  similar  displacement  field  for  the  core 
is  assumed  as  follows: 


u‘'(x,y,z,t) 


+ 


(4.3) 


6 dx 
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v^{x,y,z,t)  = 


Vq  . z 


C,c- 


dw 


\ 


0 


dy 


/ 


z^  3(|) 


z^  a 


aij;  ac 

+ 


6 dy 


y dx  dy 


/ 


w \x,y,z,t)  = Wg  + z 4)  - 


aij;  ac 

+ 


V 


dx  dy 


y 


(4.4) 


(4.5) 


where  c,  = and  Cj  = E/G^^,  Mq>  ^o>  ^o’  C and  <|)  are  functions  of  x,  y and 

t. 

The  kinematic  relations  assuming  small  linear  deformations,  take  the 
following  form: 


(4.6) 


(4.7) 


(4.8) 


where  j = t,  b correspond  to  k = 1,2  respectively. 
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The  kinematic  relation  for  the  core  is: 


y1-  = c,  ^{x,y,t) 


di]f(x,y,t)  dC,(x,yj)^ 


dx 


dy 


/ 


(4.9) 


The  compatibility  conditions  of  deformation  that  the  face-core  interface 
shear  strain  is  proportional  to  the  face-core  interface  shear  stress  and  the 
continuity  conditions  may  be  written  in  the  following: 

At  the  upper  interface  ( z = h ) 


^x-j  - k‘(u‘  - M 0 


(4.10) 


w = w ‘ 


At  the  bottom  interface  ( z = 0 ) 


'^‘xz2  " k ‘(u‘'  - M*) 


<z2-k‘(v^ -V^’) 


(4.11) 


w‘"  =w‘’ 


where  t'  ,,  '^yzi  the  shear  stresses  at  the  upper  and  lower  interface 

layers  respectively,  and  k'  = k{(,>)(\+ir],X^))  is  the  frequency-dependent  complex 
shear  modulus  of  interface  layers. 
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The  constitutive  relations  for  the  laminated  face  layers  are  derived  by 
considering  force  and  moment  resultants  at  the  top  and  bottom  faces  separately: 


(4.12) 


where  again  j = t,  b correspond  to  ^ = 1,2  respectively.  Based  on  stress-strain 
relation  of  plane  stress  state  with  respect  to  the  laminate  axes  oriented  at  an  angle 
with  respect  to  the  ply  orthotropic  axes  in  the  plane  of  the  laminate 


f > 

r 

c' 

'^ll 

^12 

'^is 

► = 

^22 

^23 

sym. 

^33 

» xy 
k -^7 

(4.13) 


the  constitutive  equations  for  the  face  laminates  will  then  be  of  the  form 


N. 


N. 


j 


./ 


J 


xy 


M 


M 


J 


J 


M 


J 


xy 


A J 4 J A n ^ n 
^11  ^12  ^16  ^11  ^ 


J 


A 

^22  ^26 


7 


sym. 


D 7 jy  J n 
^12  ^22  ^26 


7 

12 

./ 


^16  ^26  ^66 


./ 


K D 


J 


B 


J 

16 


7 


7 


7 


12  ^16 


7 


D 

^22  ^26 


D 


J 

66 


/' 


J 


y 

• xy 


k!. 


J V 


xy 


(4.14) 


where  the  A,  B,  D matrices  can  be  written  as 
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A 


j 

mn 


f 


C'  dz, 

mn  k 


h->'iP- 


B 


J 

mn 


f 


C ' z dz 

mn  k k 


(m,n  = 1,2,3) 


{m,n  = 1,2,3) 


(4.15) 


D 


J 

mn 


/ 


zldz, 

mn  k 


(m,n  = 1,2,3) 


-*  - 


and  Ky  and  are  the  curvatures  of  middle  surfaces  of  face  layers  which  can 
be  written  as 


(4.16) 


where  j = t,  b correspond  to  A:  = 1,2  respectively. 
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The  stress-strain  relations  for  the  core  are 


= K e. 


c c 


L-  ✓ — t ^ ^ 

T = G Y 

jcr  xr  • jcr 


(4.17) 


c ^ c c 

T = G Y 

vr  vz  * j:r 


Using  above  kinematic  relations  and  constitutive  equations,  we  can 
eliminate  Uq,  Vq,  Wq  and  <f)  from  the  compatibility  equations  (Eqns  (4.10)  and 
(4.11))  and  obtain  two  independent  compatibility  equations  in  the  following: 


h dWj  + h dw. 

W j ~ 1^2  ''' 


dx 


G 

2—^h 


\ 


k' 


c,  ij; 


/ 


d 


dx 


di|r  dC 


12  dx 


dx  dy 


= 0 


(4.18) 


V - V + 
1 2 


+ h dW|  h * h dw 
+ 


dy 


G 


\ 


k‘ 


h 


c.C 


/ 


dy 


d 


12  d;; 


dij;  dC 


dx  dy 


= 0 


/ 


(4.19) 


Correspondingly,  two  Lagrange  multipliers  X„  Xy  with  two  interfacial  effect 
terms  <p^,  <p^  are  introduced  into  the  total  potential  energy  equation.  The  physical 
meaning  of  Lagrange  multipliers  X^,  Xy  are  the  face-core  interface  shear  forces  in 
the  x,  y directions,  and  <p^,  <p.  can  be  regarded  as  the  strain  energy  contribution  of 


the  adhesive  layer  to  the  total  potential  energy. 
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4.2  Equations  of  Motion 


The  equations  of  motion  and  corresponding  boundary  conditions  can  be 
also  derived  by  using  the  Hamilton’s  principle.  The  total  kinetic  energy  of  the 
sandwich  plate  can  be  expressed  as  in  terms  of  eight  basic  functions, 
w,,  Wj,  i];  and  g as  follows: 


\ riJ  dw,  hf  dw, 

T ■ ^/  ]p, />,[(»,)" . (v/  • • K)^l 


12  dx  12  dy 


hi  dw  hl^  dw 


12  dx  12  dy 


+ 


(4.20) 


where  p,,  P2  and  p are  the  density  of  the  upper  face,  lower  face  and  the  core, 
respectively,  ()  denotes  a first  derivative  with  respect  to  time  and  uXx,y,z.t), 
vXx,y,z,t),  wXx,y,z,t)  are 


U = M,  + 


c G h dw 

I I ^ ^ ■ 

^ - —— — * ^ ^1  ■ 


k' 


2 dx 


dw. 

dx 


\ 


■ / 

1 

2 

[h[ 

dw^  dw^ 
dx  dx 


\ 


h d 
2 dx 


di|j  ^ dC 
dx  dx 


\ 


/ 


(4.21) 


+ 


6 dx 


/ 


\ 


dx  dy 


\ 


/ 
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c,G 

c 2 yr 

V = ^ ^ 


k‘ 


h.  dw. 

c - 

2 ay 


c,C- 


aw. 

dy 


\ 


z' 

1 

\ 

aW|  dw^ 

2 

h 

1 ay  ay  j 

h a 


2 dy 


at  ac 

+ 

ax  dx 


\ 


/ 


(4.22) 


a 


6 ay 


^ dx  ay 


\ 


/ 


W = W-  + z 

c 2 


W,  - W3 

'at  dC 

z' 

'at  ac^ 

h 2 

^ dx  dyj 

2 

h 

(4.23) 


The  total  potential  energy,  U,  consists  of  the  contributions  from  both  face 
layers  f/,  and  U2,  core  shear  strain  t/3,  core  normal  U4,  the  external  potential  term 
V due  to  applied  loading  and  the  effects  of  nonrigid  adhesive  interface  layers, 
and  <py 

f/  = t/,  . G,  . C/3  . t/,  . F . (p^  . (p^  (4.24) 

For  simplicity,  we  consider  symmetric  laminated  face  layers  where  = 0, 

„ = 0 and  = 0 (/  - I,  2;  m = n - 1,  2,  6)  in  Eqn  (4.14).  The  face 

strain  energy  contributions  then  are 
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U.  = 


if-f 

2 Jo  Jo 


/ \ 

f; 

A ' 

awj 

^ dx  ) 

+ 


12 


2^.^— 


\ 2 

av,  ^ 


ax  a>” 


dy/ 


/ 


.4 


./ 

66 


\ 


a_y  ax 


\ 


+ 


Z) 


./ 

11 


, \ 
a^w. 


/ 


ax^  / 


2D 


, a^w,  a^w, 

J k k 


12 


ax  ^ a_y  ^ 


\ 

d^w. 


dy 


^ 4d,; 


, \ 
aSv, 

k 


\ / 


dxdy/ 


dxdy 


in  which  j = t,  b when  k = 1,  2.  The  shear  strain  eontribution  of  the 


U.  = 


.G>,cy 


dx 


The  transverse  normal  strain  contribution  of  the  core  is 


U.  - ["E 

2 Jo  Jo 


(w,  -w/  /^3 

+ t + ) 

h 12  dx  dy 


dxdy 


The  contribution  of  face-core  compatibilities  are, 


= \/"/ 
./o  Jo 


+ h aW|  h + h dw 

u,  - U.+ + 

2 dx  2 dx 


G 

2—  . h 
k' 


\ 


c,  ij; 


a 


/ 


12  dx 


dx  dy 


\ 


/ J 


dxdy 


(4.25) 


core  is 
(4.26) 


(4.27) 


(4.28) 
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V - V + 
1 2 


+ h dw  h ^ h dw. 
2 dy  2 dy 


(4.29) 


G 

2^  . h 


\ 


d 


3i|;  3C 

+ 


/ 


12  a;; 


dx  dy 


/ J 


dxdy 


The  work  done  by  the  external  force  is 


dxdy 

Jo  Jo 


(4.30) 


where  q = q(x,t)  is  transverse  distributed  dynamic  loading  at  the  surface  of  top 
face  layer. 

Taking  variation  in  the  Hamilton’s  equation  with  respect  to  «i, 

Wj,  i|;,  g,  Aj.  and  A^,,  ten  equations  of  motion  and  boundary  conditions  are  obtained 
after  integration  by  parts  in  x and  t,  and  some  algebraic  manipulations  and 
substitutions.  It  is  noted  that  homogeneous  conditions  for  the  displacement  and 
the  velocity  at  ? = t,,  ?2  assumed.  The  number  of  governing  equations  is 
reduced  from  ten  to  eight  by  eliminating  the  and  A^.  The  governing  equations 
can  be  written  in  the  matrix  form: 

[dM]{u)  . [d/:]{f/}  = {F}  (4.31) 


where 


{U)  = {w^,  W2,  u^,  Uj,  V,,  V2,  aV 
{F}  = {q,  0,  0,  0,  0,  0,  0,  0}^ 
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[SM]  — the  mass  matrix  differential  operator 


— the  stiffness  matrix  differential  operator 


It  is  noted  that  [3^,  {U}  and  {U}  are  complex  valued. 

The  boundary  conditions  at  the  end  of  the  beam  can  be  obtained  in  a similar 
way  to  the  beam  theory.  Hereby  only  a type  of  simply  supported  boundary 
conditions  are  given  in  the  following, 
at  the  top  face  layer,  a:  = 0 and  x = a: 


N'  = V = M‘  = w = 0 


(4.32) 


y = 0 and  y = b: 


n;  - = m;  -w^-o 


(4.33) 


at  the  bottom  face  layer,  x = 0 and  x = a: 


(4.34) 


y = 0 and  y = b: 


v=N  = M*  = w = 0 

2 xy  y 2 


(4.35) 


at  the  core,  x = 0 and  x = a: 


= 0 (<_.  = 0) 


(4.36) 


y = 0 and  y = b: 


c = 0 (t;  = 0) 


(4.37) 
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4.3  Solution  and  Discussion 


The  solution  for  a sandwich  plate  with  simply-supported  boundary 
conditions  can  be  assumed  as 


00  00 


{U(x,i)}  . 


. mix  . muy 
w,  sin sin 

inm  1 

a b 


. mux  . muy 
sin sin ^ 

mm  7 

a b 


mix  . m'Ky 
u,  cos sin ^ 

Inm  1 

a b 


. mix  mny 
w.  sin cos ^ 


2nm 


n=0  m-0 


a 


. mix  mny 
V,  sin cos ^ 


\nm 


a 


nnx  . mny 
V.  cos sin 

inm  1 

a b 


, . nnx  . mny 

lU  sin sin 

' nm 


a 


^ . nnx  . mny 

C sin sin 

^ nm  7 

a b 


(4.38) 


where  «2„^,  v,„^,  V2„„,  i|;„„  and  are  the  amplitude  coefficients 

and  n is  the  wave  number.  This  solution  fully  satisfies  the  boundary  conditions 
of  a plate  where  faces  are  simply  supported  and  the  vertical  deflection  of  the  core 
through  the  height  is  prevented  at  four  edges.  The  constants  W2„„,  «2„m, 

v,„„,  V2„;„,  and  <;„„  are  determined  by  the  solution  of  the  eigenvalue  problem. 

As  an  example,  a sandwich  plate  consisting  of  identical  laminated  top  and 
bottom  face  layers  with  each  [0/±45/90]  lamination  and  isotropic  core  material. 
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The  following  material  properties  are  used  for  the  graphite/epoxy  lamina, 

E,  = 181  Gpa,  E2  = 10.3  Gpa,  G,2  = 7.17  Gpa 
v,2  = 0.28,  p = 1.58x10-^  N-sVm'' 

The  loss  factor  of  the  core  and  interface  layer  are  assumed  to  be  t|^  = 0.3  and  ti^ 
= 0.5,  respectively.  The  effective  modulus  of  face  layer  is  referred  to  as  E^. 

The  effect  of  interface  shear  modulus  on  plate  structural  damping  is  shown 
in  Figure  4.2.  The  similar  trend  to  the  beam  problem  can  be  observed.  For  the 
fixed  values  of  faces  and  core  moduli,  the  interface  with  some  small  value  of 
stiffness  kh/G^  results  in  a considerable  increase  in  system  damping  properties. 
This  can  be  explained  by  the  fact  that  the  smaller  value  of  k makes  E,/k  larger, 
resulting  into  larger  shear  stress  concentration  and  energy  dissipation  which  causes 
enhancement  in  damping.  Again,  a trade-off  between  damping  increase  and 
stiffness  decrease  must  be  considered. 

In  determining  the  resonance  frequency,  it  is  observed  that  the  mode 
corresponding  to  the  lowest  frequency  is  of  predominantly  flexural  type,  the 
corresponding  values  of  all  displacement  ratios  are  small.  There  are  other 
additional  families  of  modes  such  as  of  predominantly  extensional  type  and 
thickness  shear  type.  In  Table  4.1,  the  resonance  frequencies  of  the  predominantly 
flexural  mode  are  given  for  comparison  of  the  cases  with  inertia  effects  of  all 
assumed  displacements  and  with  only  transverse  inertia  effect.  The  non- 
dimensional  resonance  frequency  is  defined  as 
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khlG'^ 


Figure  4.2  Effect  of  Storage  Shear  Modulus  of  Interface  on  Loss  Factor 

of  Sandwich  Plate 
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(4.39) 


It  can  be  seen  that  the  transverse  inertia  effect  is  the  most  important  effect  on  the 
resonance  frequencies,  and  other  inertia  effects  are  negligible  for  this  mode. 


Table  4.1  Resonance  Frequencies  of  Sandwich  Plate 


Modal 

number 

Non-dimensional  frequency  QxlO^ 
for  predominantly  flexural  mode 

n 

m 

Transverse  inertia 
terms 

All  inertia 
terms 

1 

1 

1.251 

1.248 

1 

3 

11.37 

11.31 

2 

5 

69.76 

69.59 

5 

5 

145.19 

144.91 

4 

6 

232.64 

232.14 

3 

8 

294.44 

293.67 

a = b,  h/a  = 0.15,  h,/h  = h,/h  = 0.1,  G^E'  = 0.4,  = 0.01 

p/p,  = 0.2,  kh!G[  = 1000. 

Figure  4.3  shows  the  variation  of  sandwich  plate  system  loss  factor  with  the 
ratio  of  core  shear  modulus  to  face  effective  modulus.  It  is  found  that  the  system 
loss  factor  of  the  lowest  frequency  first  increases  with  core  shear  modulus  and 
then  decreases,  depending  on  the  different  aspect  ratio  of  the  plate,  /i,/a. 
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Figure  4.3  Effect  of  Core  Shear  Modulus  on  Loss  Factors  of  Sandwich  Plate 


CHAPTER  5 

VIBRATION  OF  DELAMINATED  SANDWICH  BEAM-PLATES 


In  the  Interface  model  of  sandwich  beams  discussed  in  Chapter  3,  there  are 
two  extreme  interface  conditions,  i.e.,  the  case  of  perfect  bond  between  the  face 

and  the  core  when  the  interface  adhesive  stiffness  k* — and  the  case  of 
delamination  or  interface  crack  when  k* — >-0.  Delaminations  are  usually  formed 

either  due  to  imperfections  in  the  fabrication  processes  or  under  various  external 
conditions  or  loadings  that  may  occur  in  the  operational  life  of  a composite 
structure.  Delaminations  will  affect  the  strength,  stiffness  and  integrity  of  the 
composite  structure,  and  may  thus  cause  the  change  of  structural  vibration 
response.  In  this  chapter,  a simple  beam-plate  model  is  developed  for  the 
investigation  of  the  effect  of  delaminations  on  the  vibration  behavior  of  sandwich 
composites.  The  refined  vibration  theory  of  sandwich  beams  presented  in  Chapter 
2 is  employed  for  the  analysis.  The  effects  of  delamination  size  and  position  in 
determining  vibration  characteristics  are  examined. 

5.1  A Beam-Plate  Model 

Consider  a typical  construction  of  sandwich  plates  consisting  of  two  faces, 
not  necessarily  identical,  metallic  or  laminated  symmetric  composites,  bonded  to 
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the  core  through  interface  layers.  In  modern  sandwich  structures,  plastic  foam  and 
non-metallic  honeycomb  cores  are  replacing  the  metallic  honeycomb  ones,  thus 
achieving  weight  reduction  and  flexibility  in  manufacturing  process.  In  such  cases 
the  transverse  flexibility  of  the  core  affects  the  overall  behavior,  the  stresses  and 
displacement  fields.  It  is  necessary  to  take  this  into  account  by  using  the  refined 
vibration  theory. 

A sandwich  plate  containing  an  arbitrary  located  delamination  (interface 
crack)  lying  between  the  upper  face  and  the  core  is  shown  in  Figure  5.1.  In  the 
case  of  a delamination  at  the  lower  interface  the  similar  analysis  should  be  used 
with  minor  modification.  The  plate  has  a constant  width  between  two  lateral  edges 
and  is  subjected  to  transverse  loading  q and  compressive  in-plane  force  P.  A 
beam-plate  model  assumes  that  the  delamination  extends  over  an  interval  a < x 
< b and  runs  across  the  whole  width  of  the  plate  as  depicted  in  Figure  5.1.  The 
delamination  is  assumed  to  exist  before  loading  and  not  to  grow  under  it.  It  is 
also  assumed  that  the  surfaces  of  delaminated  faces  and  core  are  free  of  shear 
stresses  (they  may  slip  one  with  respect  the  other)  but  can  accommodate  vertical 
normal  compressive  stresses  when  contact  between  them  exists. 

For  the  purpose  of  analysis,  the  delaminated  sandwich  beam-plate  has  been 
divided  into  three  spanwise  regions  as  denoted  by  circle  numbers  in  Figure  5.1, 
namely,  a delamination  region  and  two  integral  regions  one  on  other  side  of  the 
delamination  region.  The  delamination  region  is  itself  made  up  of  two  separate 
component  segment  above  and  below  the  plane  of  the  delamination,  joined  at  their 
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ends  to  the  integral  segments.  The  beam-plate  as  a whole  are  described  by 
equations  of  motion,  continuity  requirements  at  the  interface  crack  tips  and 
boundary  conditions  at  the  ends.  The  governing  equations  of  motion  and 
corresponding  boundary  conditions  for  the  regions  1 and  3 are  the  similar  one 
presented  in  Chapter  2.  For  the  delaminated  region  2,  the  governing  equations  of 
motion  and  the  continuity  requirements  are  derived  for  two  different  deformation 


modes  by  using  variational  procedures. 


Figure  5.1  Delaminated  Composite  Sandwich 
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5.2  Deformation  Mode  Modeling 

As  mentioned  in  Section  1.2.3,  there  are  two  different  modeling  of  vibration 
of  the  delaminated  laminated  beam-plate  with  respect  to  pre-buckled  state  and 
post-buckled  state,  respectively.  Nonetheless,  we  could  not  find  any  research  work 
on  vibration  analysis  of  sandwich  beam-plate  with  delaminations.  For  a 
delaminated  sandwich  beam,  we  can  model  the  delamination  region  in  two 
different  ways:  contact  or  non-contact  throughout  the  entire  length  of  the 
delaminated  region. 

In  the  case  of  a delamination  without  contact  at  the  face-core  interface  of 
delamination  region,  the  surface  are  free  of  stresses  any  kind.  The  non-contact 
assumption  implies  that  each  delaminated  segment,  face  or  core,  above  and  below 
the  plane  of  the  delamination  is  free  to  undergo  transverse  displacement  as  per  the 
requirements  of  equilibrium.  The  collaboration  between  the  two  is  a result  of  the 
continuity  requirements  at  the  ends  of  the  region.  The  delaminated  face  and  core 
are  thus  assumed  to  be  "free"  to  have  different  transverse  displacements,  as  shown 
in  Figure  5.2.  We  refer  to  such  a mode  of  deformation  as  the  free  mode.  This 
mode  may  also  be  considered  as  the  vibration  relative  to  the  post-buckling  state, 
since  there  is  an  opening  between  the  delaminated  face  and  core,  as  depicted  in 


Figure  5.2. 
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Figure  5.2  The  Free  Deformation  Pattern 


Figure  5.3  The  Constrained  Mode  (a)  Deformation  Pattern 

(b)  Contaet  Stresses  in  the  Delamination  Region 


If  the  vibration  is  considered  relative  to  the  unloaded  referential  state  or 
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pre-buckied  state  in  which  there  is  no  opening  at  the  crack  interface,  it  seems 
necessary  to  impose  the  constraint  of  continuity  condition  of  transverse 
displacement  at  the  delaminated  interface.  The  crack  surfaces  of  the  face  and  the 
core  are  assumed  to  be  free  to  slide  horizontally  with  respect  to  one  another  except 
at  their  ends,  but  remain  contact  in  the  thickness  direction.  Therefore,  the 
delamination  interface  are  subjected  to  vertical  compressive  stresses  but  are  free 
of  shear  ones,  i.e.,  the  contact  is  assumed  to  be  frictionless.  Such  a mode  of 
deformation  is  referred  to  in  this  study  as  the  constrained  mode.  In  this  model, 
the  assumption  of  no  natural  gap  between  the  delaminated  face  and  core  must  be 
maintained.  When  there  is  a gap  existing  in  the  undeformed  state,  the  problem  of 
vibration  of  delaminated  sandwich  beams  becomes,  in  general,  quite  complex  and 
non-linear. 

5.3  Analytical  Formulation 

For  the  two  models  of  deformation  modes,  the  governing  equations  of 
motion  and  the  boundary  conditions  for  the  delaminated  regions  are  derived  using 
the  Hamilton  principle,  similar  to  the  procedure  described  in  the  Previous  chapters. 
Our  primary  concern  is  the  delamination  with  contact,  i.e.,  the  constrained  mode 
model  because  this  is  a realistic  model,  whereas  the  limitation  of  the  free  mode 
model  is  discussed  via  some  typical  results. 
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5.3.1  The  Free  Mode  Model 


In  the  case  of  delamination  without  contact,  the  stress  free  conditions  at  the 
delaminated  interface,  i.e.,  the  upper  interface,  requires 


o,(x,  z =h,  t)  = 0 
zl(x,  z = h,t)  = 0 


(5.1) 


Recall  the  kinematic  relations  for  the  core,  Eqn  (2.6)  in  Chapter  2,  resulting  the 
following  equations 


4)(jc,t)  = 0 


4r(x,t)  = 0 


(5.2) 


The  compatibility  conditions  of  deformation  at  the  bounded  face-core 
interface,  i.e.,  the  lower  interface,  yield  the  same  equations  as  in  Chapter  2, 
at  z = 0 


u ‘^  = u' 


(5.3) 


which  give  the  relations  between  the  displacement  functions  as  follows 


wp:,t)  = w^{x,t) 


Mp(jC,/)  = U^(x,t)  - 


/?2  dW2(x,0 

2 dx 


(5.4) 
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Now  the  independent  displacement  functions  are  u^(x,t),  w^(x,t),  U2(x,t)  and 
W2(x,t).  Using  the  Hamilton  principle  again,  the  governing  equations  can  be 
derived  via  variational  calculus.  The  total  kinetic  energy  of  the  sandwich  segment 
at  the  delamination  region  can  be  expressed  in  terms  of  four  independent  functions 
as  follows 


T 


1 fM  u 


(w, )' 


+ 


12 


/ \ 
dw 


1 


\ 


dx 


+ 


f 

1 

dw- 

2 

{Uif  . 

^ 12 

2 

dx  J 

. (Wif 

+ 


p h 


(«2)' 


- (/z  + 


dw^ 

dx 


h 


2 hh^  h2 
— + + 


\ 


\ 

/ \ 

9 

> 

40 

2 

> 

/ 

dx  J 

J 

dx 


(5.5) 


The  total  potential  energy,  U,  consisting  of  the  contribution  from  both  faces,  and 
the  external  potential  terms,  can  be  written  as 


U = 


h. 


du 


, h, 


1 


El  hi 


du 


12  dx^ 


hi  di^w 

^^2  , _i_( If 


12 


(5.6) 


qw^  - (P  ‘u  ' - P ‘’u  *)6^x  -x^  \dx 


where  q = q{x,t)  is  transverse  distributed  dynamic  loading  at  the  surface  of  the 
top  face  layer,  and  P,  P are  longitudinal  compressive  forces  applied  at  upper  and 


lower  faces,  respectively. 
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Taking  variation  in  the  Hamilton’s  equation  with  respect  to  w^,  W2,  «i  and 
four  equations  of  motion  and  consistent  boundary  conditions  are  obtained  after 
integration  by  parts  in  x and  t,  and  some  algebraic  manipulations  and  substitutions. 
It  is  noted  that  homogeneous  conditions  for  the  displacement  and  the  velocity  at 
t = t,,  ?2  are  assumed.  The  four  governing  equations  can  be  written  in  the  matrix 
form: 

[dM]{u]  . [ax]{t/}  = {F}  (5-7) 

where 

{U}  = {w,,  W2,  M,,  U2V 
{F}  = {q,  0,  0,  0}^ 

[dM]  — the  mass  matrix  differential  operators 
[0^  — the  stiffness  matrix  differential  operators 

The  continuity  requirements,  at  the  end  of  the  delamination  region,  x=a  or 
x = b,  (see  Figure  5.1),  consist  of  six  geometrical  conditions  and  six  natural 
conditions. 

For  the  upper  face  layer,  we  have 

u\\x=a,t)  = u^{x=a,t) 
u^^\x  = b,t)  = u = b,  t) 


(5.8) 
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w\\x=a,i)  = 

(2) 

w^  (x=a,t) 

wf\x=b,t)  = 

wWx^b,f) 

dw^^\x=a,t) 

dx 

dx 

dw^^\x  = b,t) 

u> 

II 

dx 

dx 

t(2) 

(x=a,t)  = 

(x = a,  t) 

A2) 

,(3) 

Nj  (x-b,t)  - 

N‘  {x-b,t) 

f{2) 

{x  = a,t)  = 

(x=a,t) 

A2) 

/3) 

(x  = b,t)  = 

(x=b,t) 

A2) 

dM^  (x = a,  t) 

dM^  (x=a,t) 

dx 

dx 

A2) 

,(3) 

dM'  {x-b,t) 

dM'  (x=b,t) 

dx 

dx 

(5.9) 


(5.10) 


(5.11) 


(5.12) 


(5.13) 


where  the  superscripts  (1),  (2)  and  (3)  represent  three  regions.  It  is  noted  that  in 
the  shear  force  balance  conditions  at  the  delamination  front,  Eqn  (5.13),  the  inertia 
effects  are  neglected. 
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For  the  lower  face  layer,  we  have  the  similar  continuity  conditions, 


u\\x  = a,t)  = u?\x  = a,t) 


u^\x  = b,t)  = uf\x  = b,t) 


WjXx^a,!)  = W2\x=a,t) 


w!^\x  = b,t)  = w^\x  = b,t) 


(5.14) 


(5.15) 


= a,  t)  dw^\x = a,  t) 


dx 

dx 

dw2^\x  = b,t) 

dw2^\x  = b,f) 

dx 

dx 

b<^K 

{x  = a,  t) 

= (x  = a,t) 

K (x-b,t) 

= (x  = b,t) 

/,(2) 

{x  = a,  t) 

= (x  = a,  t) 

(x  = b,t) 

= (x  = b,t) 

dM^  (x = a,  t) 

dM^  (x = a,  t) 

dx 

dx 

dM^  (x  = b,t) 

aM;  (x=b,t) 

dx 

dx 

(5.16) 


(5.17) 


(5.18) 


(5.19) 
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where  N(  and  Ml  (j  = t,  b)  are  the  axial  forces  and  the  bending  moments  in  the 
faces,  respectively,  and  the  superscripts  (/),  / = 1,  2,  3,  denote  the  three  spanwise 
regions  as  described  in  Section  5.1.  The  constitutive  equations  relating  N{  and  M{ 
to  the  corresponding  displacements  are  given  in  Eqn  (2.16). 

There  is  only  one  continuity  condition  imposed  on  the  core  at  each  end  of 
delamination  region. 


c<'>  c® 

{x=a,i)  = (x=a,t)  = 0 

^(2)  „(3) 

(x  = b,t)  = (x=b,t)  = 0 


(5.20) 


If  the  delamination  region  starts  or  ends  at  one  of  the  beam  edges,  for  example, 
at  j:  = 0,  the  continuity  conditions  at  one  end  are  replaced  with  boundary 
conditions  which  yield, 
at  the  upper  face  layer  { x = 0 ): 


N‘  = -P‘  or  w,  = 0 


(5.21) 


(5.22) 


dj^  ^ P,/?'  aw,  p/;2 

dx  12  dx  2 

or  w = 0 


(5.23) 
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at  the  bottom  face  layer  ( x = 0 ): 


(5.24) 


(5.25) 


0X 


P2^2  ^2 
12  3x 


% 2h 

3 40  0x  15  dx  ^ 


or  ^2=0 


and  at  the  core  ( x = 0 ): 


(5.26) 


i|;  = 0 


(5.27) 


5.3.2  The  Constrained  Mode  Model 

Consider  a delamination  region  with  contact,  the  crack  surfaces  remaining 
in  contact  vertically  but  free  of  shear  stresses.  Thus,  the  shear  stresses  in  the  core 
are  null,  which  can  be  represented  as, 

il;(x, /)  = 0 (5.28) 

The  formulation  procedure  is  identical  with  that  of  the  free  mode  model,  which 
leads  to  similar  governing  equations  of  motion  and  similar  continuity  and  boundary 
conditions.  Some  differences,  however,  exist  due  to  the  contact  normal  stresses 
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at  the  delaminated  surfaces.  Such  differences  are  represented  by  the  compatibility 
condition  between  the  core  and  the  faces  at  the  upper  interface  in  the  following: 

w ‘{x,z^=--^,t)  = w ‘'(x,  z = h,t)  (5.29) 

The  compatibility  conditions  between  the  core  and  the  faces  at  the  lower  bonded 
interface  are  the  same  as  that  of  the  region  without  contact  in  Eqn  (5.3). 
Substituting  Eqns  (5.28),  (5.29)  and  (5.3)  into  displacement  fields,  Eqns  (2.1), 
(2.2)  and  (2.3),  the  following  equations  are  given. 


u 


0 


2 dx 


w 


i|;  = 0 


(5.30) 


which  also  results  in  four  independent  displacement  functions,  u^{x,t),  w,(x,f), 
U2(x,t)  and  \V2(x,t),  as  seen  in  the  free  mode  model. 

The  total  kinetic  energy  of  the  sandwich  segment  at  the  delamination  region, 
having  different  core  contribution  than  the  free  mode  model,  can  be  expressed 
in  the  form 

T = ^ ^ (5.31) 

where  the  kinetic  energy  contributions  from  two  face  layers,  T,  and  Tj,  are  the 

same  as  that  of  the  free  mode  model,  and  can  be  written  as. 


105 


(5.32) 


It  can  be  seen  in  Eqn  (5.6)  that  the  total  potential  energy  U does  not  have  the 
contribution  from  the  core  in  the  case  of  delamination  without  contact.  Since 
there  are  contact  normal  stresses  in  the  constrained  model,  a potential  energy 
contribution  from  the  core  normal  stresses  must  be  included  in  the  total  potential 
energy  U and  it  gives. 
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(5.33) 


The  four  governing  equations  of  motion  can  be  derived  in  a similar  way  as 
used  in  the  free  mode  model.  The  continuity  conditions  at  the  ends  of  the 
delamination  region  will  be  the  same  as  Eqns  (5.8)  through  (5.20).  In  the  case  of 
the  delamination  region  starting  at  the  beam  edge  {x  = 0),  the  continuity  conditions 
at  that  end  become  boundary  conditions  which  are  the  same  as  Eqn  (5.21)  through 
(5.27).  Similarly,  the  boundary  conditions  can  be  obtained  if  the  delamination 
region  ends  at  the  beam  edge  {x  = /). 
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5.4  Example  Problems  and  Discussion 

A few  numerical  example  problems  are  considered  for  illustrative  purpose. 
To  verify  the  theory  and  computer  program,  the  first  three  natural  frequencies  of 
two  different  models  were  calculated  for  a delamination  length  of  0.0001/  and 
showed  negligible  discrepancies  from  those  of  a perfect  sandwich  beam  considered 
in  Chapter  2. 

5.4.1  The  Free  Mode  Model 

The  problem  of  the  free  vibration  {q  = 0)  with  simply  supported  boundary 
conditions  can  be  solved  in  the  similar  way  discussed  in  Chapter  2 and  3 by 
assuming  a harmonic  motion.  The  natural  frequencies  and  mode  shapes  are 
calculated  for  a sandwich  beam  with  a delamination  located  symmetrically  along 
X direction.  Figure  5.4  (a)-(c)  show  the  normalized  mode  shapes  for  the  case, 
d/l  - 0.3,  hy/h  = 0.1  and  h2/h  = 0.02,  h/l  = 0.05 
E,  = Ej,  E,/E,  = 0.05,  E,/G,  = 2.5 

where  d = b - a \s  the  delamination  length.  The  transverse  displacements  W;  (i  = 
1,2)  and  coordinates  x,  z are  non-dimensionalized  as 
vv,  = w,//z,  (i  = 1,2),  X = x/l,  z = zfih+h^+h^) 

The  effect  of  the  delamination  length  on  the  natural  frequencies  has  been 
examined  for  this  case  of  example  problems.  The  first  three  natural  frequencies 
of  the  bending  mode  are  given  in  Table  5.1.  In  this  table,  the  non-dimensionalized 
frequency  is  defined  as 


X 

(a)  First  Mode 


X 


(b)  Second  Mode 
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Figure  5.4  Mode  Shapes  of  Vibration  of  a Sandwich  Beam-Plate  Based  on 

the  Free  Mode  Model, h^/h  = 0.1, h^lh  = 0.02 


Table  5.1  Natural  Frequency  Q,  of  Bending  Modes 


d/l 

0.10 

0.20 

0.30 

0.40 

0.50 

0.70 

Q, 

3.813 

3.804 

3.781 

3.755 

3.717 

3.627 

Q2 

7.652 

7.324 

6.701 

6.140 

5.704 

4.651 

Q3 

14.41 

13.64 

12.25 

10.86 

9.78 

8.26 

For  a perfect  beam,  Q,  = 3.906,  Q2  = 7.811,  Q3  = 14.84 
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(5.34) 


where  the  subscript  of  Q denotes  the  mode  number.  It  is  found  that  delamination 
has  significant  effects  on  the  natural  frequencies.  All  the  three  frequencies  are 
lower  than  those  of  the  perfect  sandwich  beam,  and  the  difference  between  the 
second  and  third  natural  frequencies  is  larger  than  the  difference  between  the  first 
and  second  ones. 

It  can  be  seen  that  in  the  case  of  the  first  two  symmetric  modes  as  shown 
in  Figure  5.4(a)  and  (b),  the  mode  shapes  of  the  face  and  core  in  the  delamination 
region  are  different  and  that  the  transverse  displacement  of  the  face  is  greater  than 
that  of  the  core  at  all  corresponding  axial  positions.  Such  a free  displacement 
mode  would  give  rise,  in  one  half  of  the  cycle  of  motion,  to  an  overlap  between 
the  delaminated  face  and  core.  This  is  a downward  displacement  of  the  face 
greater  than  that  of  the  lower  core.  In  the  case  of  the  first  antisymmetric  mode 
as  shown  in  Figure  5.4(c),  the  mode  shapes  would  have  an  overlap  over  a half  of 
the  delamination  region  during  the  complete  cycle  of  motion.  It  can  be  shown  that 
similar  overlaps  will  appear  at  all  higher  modes  in  this  free  mode  model. 
Therefore,  the  deformation  modes  given  by  the  free  mode  model  are  not  physically 
feasible,  although  mathematically  admissible  as  equilibrium  solutions. 

As  pointed  out  earlier,  the  free  mode  model  may  be  considered  as  the 
vibration  relative  to  the  post-buckling  state.  In  this  case,  the  governing  equations 
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of  motion  must  therefore  account  for  the  geometric  stiffness  produced  by  the  axial 
loads.  As  these  loads  are  themselves  functions  of  the  transverse  displacements,  the 
terms  representing  the  geometric  stiffness  in  the  governing  equations  of  motion  are 
non-linear.  These  non-linear  terms,  which  are  ignored  in  the  analysis  of  free  mode 
model,  may  become  significant  even  within  the  regime  of  small  deformation 
analysis,  depending  on  the  size  and  location  of  the  delamination.  As  a result,  the 
governing  equations  of  motion  are  valid  only  over  a very  small  portion  of  the  free 
mode  deformation  cycle  during  which  it  is  physically  feasible. 

5.4.2  The  Constrained  Mode  Model 

The  above  discussions  showed  that  the  free  mode  model  gives  vibration 
modes  with  overlaps  of  deformation  which  violate  geometric  compatibility.  In 
reality,  however,  the  tendency  of  the  delaminated  face  layers  to  overlap  on  the 
core  will  be  resisted  by  the  development  of  a contact  stress  distribution  between 
the  face-core  interface.  Such  a contact  stress  distribution  could  constrain  the 
transverse  deformation  of  these  face  and  core  to  be  identical  and  thus  ensure  the 
compatibility.  This  physically  admissible  mathematical  model  has  been  widely 
used  and  verified  by  experiments  for  vibrations  of  laminated  composites  with 
delaminations  [76,  83,  84],  where  the  identical  transverse  displacements  of  the 
delaminated  layers  were  assumed.  In  the  constrained  mode  model  of  sandwich 
beams,  the  transverse  displacements  of  upper  and  lower  faces  may  be  different, 
although  the  compatibility  condition  at  the  delaminated  interface  remains. 
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In  Figure  5.5(a)  and  (b),  the  normalized  mode  shapes  of  the  first  two 
bending  modes  are  shown  for  a simply  supported  sandwich  beam  with  the  spanwise 
location  of  the  delamination  located  at  the  midpoint.  The  system  parameters  are 
similar  to  the  free  mode  model  as  defined  in  the  following: 

h^/h  = 0.02,  h/l  = 0.05,  E,  = Ej,  E,/E,  = 0.05,  E,/G,  = 2.5 
The  figure  shows  that  the  transverse  displacements  w,  and  W2  in  the  delamination 
region,  comparing  to  the  free  mode  model,  are  very  close.  It  is  alos  noted  that  the 
second  mode  in  the  free  mode  model  disappears  in  this  case. 

Figure  5.6  shows  the  effect  of  delamination  size  on  the  non-dimensional 
natural  frequencies  of  the  first  three  bending  modes.  It  can  be  observed  that  the 
reduction  in  the  frequency  of  the  first  mode  is  negligible  for  small  delamination 
size  up  to  70%.  The  weakening  effect  of  the  delamination  on  the  second  mode 
and  the  third  mode  is  much  larger  than  that  on  the  first  mode  frequency.  There 
is  a substantial  drop  in  the  frequency  values  even  for  the  small  delamination  on 
the  high  modes. 

Figures  5.7  and  5.8  show  the  effect  of  spanwise  location  of  the  delamination 
on  the  frequencies  for  the  first  two  bending  modes  of  sandwich  beams  with  simply 
supported  boundary  conditions,  respectively.  It  can  be  seen  that  the  spanwise 
location  of  the  delaminations  has  strong  effect  on  the  frequency  values.  Such 
effects  may  be  explained  by  the  fact  that  the  delamination  moves  to  regions  of 
high  shear  force  or  regions  of  high  curvature.  In  some  specific  locations  the  effect 
of  delamination  on  the  frequencies  may  become  a minimum. 
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X 

(a)  First  Mode 


Figure  5.5  Mode  Shapes  of  Vibration  of  a Sandwich  Beam-Plate  Based  on 

the  Constrained  Mode  Model,  d/l  = 0.3,  h^/h  = 0.02,  h/l  = 0.05 
E,  = Ej,  E^E,  = 0.05,  E,/G,  = 2.5 
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Figure  5.6  Non-dimensional  Frequency  versus  Delamination  Size  for  the 

First  Three  Bending  Modes 
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Figure  5.7  Effect  of  Spanwise  Location  on  the  Fundamental  Frequency 
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Figure  5.8  Effect  of  Delamination  Spanwise  Location  on  the  Frequency 

of  the  Second  Mode 


CHAPTER  6 

VIBRATION  OF  LAMINATED  BEAMS  WITH  DELAMINATIONS 


Free  vibration  of  laminated  plates  with  delaminations  has  been  examined  by 
many  researchers  analytically  and  experimentally  [73-84].  A more  general  theory 
for  modeling  of  delamination  in  composite  laminates  is  based  on  a layer-wise 
displacement  field.  The  finite  element  model  of  layer-wise  theory  is  not  only 
providing  an  adequate  framework  for  the  delamination  analysis  of  laminated 
composites,  but  also  capable  of  dealing  with  complex  geometries  of  delaminations 
(e.g.,  multiple  delaminations)  and  structures.  However,  the  computational  cost  of 
the  analysis  is  very  expensive,  and  thus  makes  it  unattractive  for  delamination 
analysis  in  engineering  practice  when  the  number  of  layer  is  large.  We  try  to 
implement  this  layer-wise  finite  element  model  in  vibration  analysis  with  emphasis 
on  developing  parallel  algorithms  of  finite  element  analysis  to  minimize  both  the 
computational  cost  and  communication  cost. 

6.1  The  Laver-Wise  Finite  Element  Model 

In  delamination  analysis  of  laminated  composites,  the  layer-wise  theory 
provides  a compromise  between  the  3-D  theory  and  conventional  plate  theories. 
Conventional  theories  based  on  a single  continuous  and  smooth  displacement  field 
through  the  thickness  of  a composite  laminate  give  poor  estimation  of  the 
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interlaminar  stresses.  With  layer-wise  continuous  representation  of  displacements 
through  the  thickness,  it  predicts  the  interlaminar  stresses  very  accurately. 
Moreover,  it  also  has  the  advantage  of  all  plate  theories  in  the  sense  that  it  is  a 
two-dimensional  theory,  and  doesn’t  suffer  from  aspect  ratio  limitations  associated 
with  three-dimensional  finite  element  models. 

For  modeling  delamination  of  laminated  composites,  the  layer-wise  theory 
can  be  extended  to  model  the  kinematics  of  a layered  plate  with  provision  of 
delamination.  It  is  implied  that  the  same  displacement  distributions  in  individual 
layer  is  capable  of  representing  displacement  discontinuity  condition  at  interfaces 
between  layers.  The  present  study  on  vibration  behavior  uses  the  same  approach 
presented  earlier  by  Barbero  and  Reddy  in  the  buckling  analysis  [67].  According 
to  Barbero  and  Reddy,  the  in-plane  displacements  are  assumed  to  have  linear 
through-the-thickness  variation  within  each  ply,  i.e.,  the  global  approximation  of 
the  form  of  in-plane  displacements  is  refined  to  a layer-by-layer  approximation. 
High  order  shear  deformable  theories  that  allow  for  quadratic  or  higher  order 
variations  of  the  in-plane  displacements  are  possible  within  the  context  of  the 
layer-wise  theory,  but  the  layer-wise  first  order  shear  deformable  assumption  is 
sufficient  for  the  present  study.  For  simplicity,  it  is  assumed  that  the  delamination 
extends  through  the  entire  width  and  there  are  no  variations  in  the  y direction,  as 
discussed  in  the  beam-plate  model.  Furthermore,  the  location  and  length  of  the 
delaminations  are  assumed  to  be  known  and  fixed.  As  a result,  the  beam  element 
of  the  layer-wise  finite  element  model  is  employed. 
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Based  on  these  assumptions,  the  global  layer-wise  displacement  terms, 
u{x,z,t)  and  w{x,z,t),  which  include  delaminations  can  be  represented  in  the 
following  form 

N D 

u(x,z,t)  = u^(x,t)  + /)  + ^2  H '{z)U \x,t) 

j-\  /.I 

(6.1) 

D 

w(x,z,t)  = w^(x,t)  t H\z)W'{x,t) 

/.  I 

where  the  step  functions  H are  computed  in  terms  of  the  Heaviside  step  functions  H 
as: 


H‘(z) 

- H{z-z)  - 1 

for  z ^ z. 

l 

H\z) 

= H(z  - z)  = 0 

for  z<z. 

i 

(6.2) 


In  Eqn  (6.1),  the  first  terms,  Uq  and  Wq,  are  the  displacement  functions  of  a 
reference  surface;  u!  is  the  layer-wise  in-plane  displacement  function  in  the  x 
direction;  and  4>'  is  a local  linear  Lagrangian  interpolation  function.  Displacement 
functions  that  represent  relative  in-plane  slippage  and  relative  transverse  opening 
displacements  across  the  ith  delaminated  interface  are  represented  by  IP  and  W, 
respectively.  Using  the  step  functions  H\  we  can  model  any  number  of 
delaminations  through  the  thickness;  the  number  of  D is  equal  to  the  number  of 
delaminations  considered.  Each  ply  is  assumed  to  be  inextensible  in  the  z 
direction  and  the  in-plane  slippage  is  assumed  to  occur  across  frictionless 
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Figure  6.1  Laminated  Composite  with  a D elamination 


interfaces.  The  reference  surface  is  taken  as  the  bottom  of  the  laminate,  as  shown 
in  Figure  6.1. 

Following  the  small  deformation  theory,  the  strain-displacement  relations 
and  constitutive  equations  can  be  obtained  for  the  layer-wise  displacement  field. 
By  taking  the  first  variation  of  the  total  kinetic  energy  and  potential  energy  or 
using  the  principle  of  virtual  work,  we  can  arrive  at  the  governing  equations  of 
motion.  To  develop  the  corresponding  finite  element  model,  the  generalized 
displacements  (Uq,  Wq,  uf,  U',  W)  are  expressed  over  each  element  as  a linear 
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combination  of  the  interpolation  function  ij;,  and  the  nodal  values  (Mq*,  vVq*,  u{,  U[, 
W[)  as  follows, 

(6.3) 

k = 1 

where  n is  the  number  of  nodes  per  element.  Then,  substituting  this  model  into 
the  governing  equations  of  motion,  we  obtain  the  following  finite  element 
equations: 
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where 
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{A}^  = {«oi,  Woi,  Uq„,  Wq„} 


{A-'y  = {u{,  u{,  ui) 


{A'}^=  {U\,  w;, U‘„  w;,} 


and  the  submatrices  [m“],  [wj'],  [mil]  and  [A:"],  [ky], 

[ky],  [ky^],  [ky^],  [kyy,  [kH]  with  i,  j = 1,  and  r,  s = \,  D are  the  mass 
matrices  and  stiffness  matrices.  There  are  (2+A^+2D)  displacement  variables  in 
each  node  of  the  one  dimensional  finite  elements  used  in  representing  the  2-D 
geometry.  The  load  vectors  {q},  {^},  and  {q'y  are  similar  to  the  forms  of  {A}, 
{A^}  and  {A^}.  In  the  case  of  free  vibration,  all  load  vectors  are  zero  and  the 
following  eigenvalue  problem  can  be  formulated  when  the  harmonic  motion  is 
assumed: 


(-«=[m].[a:])|a).  {0) 


(6.5) 


6.2  Parallel  Algorithm  Design  for  Finite  Element  Analysis 

As  stated  earlier,  one  of  the  disadvantages  of  layer-wise  finite  element 
modeling  is  the  expensive  computational  cost,  particularly  with  large  number  of 
layers  and/or  multiple  delaminations.  The  computational  power  of  the  new  parallel 
processing  computer  system  offers  an  oppotunity  to  reduce  computational  times 
for  analysis  so  that  the  finite  element  analysis  can  be  more  fruitfully  used  by  the 


designers. 
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It  is  now  clear  that  the  analysis  and  design  of  practical  engineering  systems 
in  the  future  will  make  extensive  use  of  concurrent  processing  on  computers  with 
multiple  processors.  Parallel  processor  computers  consisting  of  multiple 
processors  linked  to  one  another  and  running  subtasks  concurrently  are  becoming 
widely  available  for  large-scale  scientific  and  engineering  computations.  These 
computers  are  being  increasingly  used  to  provide  the  reduction  in  effective 
calculation  time  (or  elapsed  time)  needed  to  carry  out  many  real-time  applications. 
A new  area  of  research  in  computational  mechanics  is  based  upon  the  unique 
architectural  features  of  parallel  processing  computers.  However,  there  is  still 
a big  gap  between  well-developed  parallel  algorithm  analysis  and  their  application 
in  the  advanced  engineering  programs.  Comparing  to  the  implementation  of 
sequential  code,  parallel  algorithm  design  in  the  computational  structural 
mechanics  requires  a thorough  understanding  of  the  nature  of  the  problem  and  a 
considerable  investment  in  analyzing  the  data  flow  in  the  solution  procedure  itself 
as  well  as  the  knowledge  of  parallel  algorithm  theory. 

In  this  study,  we  try  to  develop  approximately  optimized  parallel  algorithms 
in  the  layer-wise  finite  element  analysis  of  delamination  problems.  Although  only 
such  a finite  element  model  has  been  implemented,  it  is  expected  that  our  approach 
has  the  potential  of  handling  various  numerical  schemes  for  different  purposes  of 
finite  element  analysis,  and  to  be  as  much  as  possible  hardware  independent.  It 
is  noted  that  efficient  non-numerical  algorithms  such  as  the  concurrent  assembly 
of  a mass  matrix  or  a stiffness  matrix  and  fast  schemes  for  the  solution  of  a set  of 


123 


linear  equations  may  be  related  with  specific  parallel  computers  or  architectural 
issues. 

A number  of  special  strategies  can  be  used  to  increase  the  degree  of 
parallelism  in  finite  element  computations.  These  strategies  are  applications  of 
"divide  and  conquer",  based  on  breaking  a large  (and/or  complex)  problem  into  a 
number  of  smaller  (and/or  simpler)  subproblems  which  may  be  solved 
independently  on  distinct  processors.  The  degree  of  independence  of  the 
subproblems  is  a measure  of  the  effectiveness  of  the  algorithm  since  it  determines 
the  amount  and  frequency  of  communication  and  synchronization. 

Two  parallel  algorithms  in  finite  element  analysis  are  discussed  in  the 
following:  domain  decomposition  and  optimal  mapping. 

6.2.1  Domain  Decomposition 

The  basic  idea  of  domain  (or  spatial)  decomposition  is  to  divide  the  domain 
into  a number  of  (possibly  overlapping)  regions.  The  initial/boundary-value 
problem  is  decomposed  into  one  that  involves  solution  of  initial/boundary-value 
problems  on  subproblems,  thereby  introducing  spatial  parallelism.  Given  a number 
of  available  processors  N , it  is  attractive  to  decompose  a finite  element  domain 
R into  the  same  number  of  subdomains,  because  this  allows  individual  subdomain 
operations  to  be  performed  concurrently  on  separate  CPU’s.  In  this  way,  we  are 
decentralizing  the  organization  of  the  finite  element  software  system  in  the  sense 
of  data-flow.  Each  processor  pj  is  mapped  onto  a subdomain  R:  in  which  the  basic 
geometrical  and  physical  properties  of  those  elements  and  some  additional  control 
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parameters  are  stored.  In  stead  of  seeking  a perfect  match  between  the 
communication  structure  of  the  hardware  and  the  communication  requirements  of 
subsequent  algorithms,  we  develop  an  approximate  optimal  algorithm  with 
generality  while  we  try  to  maintain  a balanced  utilization  of  a given  hardware  and 
the  software  architecture  to  be  designed. 

We  propose  a hardware-independent  domain  decomposition  algorithm  to 
fulfill  the  requirements  of  generality  in  arbitrary  discretization,  balanced 
computational  loads  and  minimal  communication  and/or  synchronization. 
Therefore,  it  must  yield  a set  of  balanced  subdomains  to  make  the  overall 
computational  load  as  evenly  distributed  as  possible  among  the  processors,  and  it 
must  minimize  the  number  of  interface  nodes  to  reduce  the  communication  cost. 

We  define  an  internal  boundary  element  within  a subdomain  Rj  as  an 
element  that  is  connected  to  an  internal  boundary  node,  which  is  the  same  as  an 
interface  node  of  the  global  domain  R.  The  domain  decomposition  algorithm  will 
provide  the  local  information  on  a list  of  elements,  internal  boundary  nodes  and 
internal  boundary  elements  for  each  subdomain  Rj,  and  global  information  on  a list 
of  the  subdomains  that  connect  to  each  interface  node,  the  connectivity  table 
(adjacency  graph)  for  the  subdomain  Rj  and  the  colored  graph  associated  with 
internal  boundary  elements.  For  a distributed  multiprocessor  computer,  as  an 
example,  only  the  node  geometry  and  element  properties  of  subdomain  R:  need  be 
stored  in  the  RAM  of  processor  pj  if  pj  is  assigned  to  Rj.  The  algorithm  then 
makes  use  of  a list  of  elements  in  each  subdomain  to  load  initially  each  processor 
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with  the  appropriate  data.  The  communication  structure  of  the  problem  is 
established  by  the  subdomain  connectivity  table.  Two  subdomains  are  considered 
as  connected  if  they  share  at  least  one  interface  node.  Then,  the  list  of  the 
subdomains  that  connect  to  each  interface  node,  the  connectivity  table  and  the  list 
of  internal  boundary  nodes  determine  the  number  and  content  of  message  to  send 
and/or  receive  for  each  processor. 

The  key  here  is  to  properly  deal  with  the  internal  boundary  elements 
because  they  constitute  a source  of  inherent  serializing  computations  and  potential 
memory  conflict  during  parallel  computations  that  are  carried  on  at  the  element 
level.  It  can  be  seen  that  each  of  these  elements  in  a subdomain  Rj  has  a 
corresponding  element  in  some  subdomain  R:  connected  to  it  through  a set  of 
interface  nodes.  We  offer  a possible  solution  to  minimize  the  serializations  in 
these  critical  regions.  We  eliminate  internal  boundary  elements,  based  on  the  local 
information,  from  their  original  subdomains  and  treated  as  a set  of  elements 
defining  a single  subdomain  R*.  Then,  a usual  coloring  technique  is  applied  to 
decompose  D*  into  a set  of  subdomains  D*  such  that  each  element  corresponds  to 
a single  color,  i.e.,  no  two  elements  are  connected.  For  each  color,  because  they 
are  mutually  disjointed,  the  elements  may  be  assigned  arbitrarily  to  the  Np 
processors.  This  algorithm  has  been  implemented  in  the  vibration  analysis  of 
delaminated  sandwich  beams  with  a regular  finite  element  mesh.  In  our  practice, 
this  domain  decomposition  algorithm  can  improve  the  computational  time  at  the 
range  of  30%  ~ 40%  in  the  case  of  using  16-processors. 
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6.2.2  Optimal  Mapping  of  Irregular  Finite  Element  Domains 

For  the  case  of  irregular  finite  element  domain  problems  such  as  structures 
with  delamination  and  fracture  or  complex  geometry,  the  good  strategy  of  mapping 
solution  domain  into  parallel  processors  is  no  longer  obvious.  An  arbitrary 
clustering  of  elements  to  form  subdomains  equal  to  the  number  of  available 
processors  will  undoubtedly  result  in  imbalanced  workload  distribution  among 
processors.  The  resources  required  by  each  processor  depend  not  only  on  the 
computation  cost  associated  with  the  elements  falling  within  its  subdomain,  but 
also  on  the  cost  of  transferring  the  subdomain  boundary  information  to  other 
processors  containing  contiguous  subdomains.  And  even  if  an  equal  number  of 
elements  are  allocated  per  processor,  excessive  communication  costs  may  strongly 
outweigh  any  parallel  processing  gains.  We  use  another  approximate  optimal 
mapping  strategy  which  guides  the  process  of  decomposing  complex  solution 
domains  into  appropriate  subdomains  so  as  to  balance  the  workload  required  by 
each  processor. 

First,  the  cost  function  C as  our  criterion  for  selecting  as  optimal  mapping 
is  defined  in  the  following 

c = E(w,f  (6.6) 

^ ^ comp  P’  ^ 

where 

w,  = computational  workload  associated  with  elements  residing  in  subdomain  i. 


assigned  to  processor  i. 
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Dpq  (Dij)  = communication  distance  between  two  finite  elements  p and  q 

residing  in  two  different  processors  i and  j. 
tcomm’  Icomp  ^ typical  communication,  computation  time  per  element,  respectively. 

The  first  term  of  the  cost  function  serves  to  equalize  the  computational 
workload  of  various  processors.  The  second  term  in  the  cost  function  represents 
the  communication  cost,  with  the  ratio  (tcomm/t^onip)  introduced  to  weigh  the  relative 
amount  of  communication  to  computation  required  per  element.  This  ratio 
depends  not  only  on  the  nature  of  the  computational  problem  being  modeled,  but 
also  on  the  relative  speed  of  the  processors  in  communicating  one  floating  point 
number  from  one  processor  to  its  nearest  neighbor  to  its  speed  in  executing  one 
floating  point  operation  on  that  number.  As  such,  (tcomm/tcomp)  is  both  problem 
dependent  and  machine  dependent. 

Clearly,  the  cost  function  C will  assume  different  values  for  different 
decompositions  of  the  solution  domain  into  subdomains.  And  even  if  one  assumes 
that  the  amount  of  computations  associated  with  each  finite  element  is  constant 
regardless  of  the  nature  of  domain  decomposition  and  processor  assignment,  it  is 
important  to  note  that  the  interelement  communication  costs  will  vary.  Thus,  the 
total  cost  of  processing  each  element,  and  therefore  all  elements,  will  depend  on 
the  relative  allocation  of  the  n-elements  to  the  A-processors. 

It  is  theoretically  possible  through  enumeration  to  exactly  determine  the 
best  assignment  policy  which  will  minimize  C.  However,  to  do  this  would  require 
one  to  evaluate  and  sort  C as  many  times  as  the  number  of  possible  ways 
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(combinations)  S for  clustering  n-elements  into  N-subdomains.  The  number  of 
these  combinations  S can  be  found  from  the  theory  of  the  classical  "Occupancy" 
problem.  Unfortunately,  this  method  is  practical  only  for  the  simple  cases  when 
n and  N are  small  (<10). 

Instead,  we  use  suboptimal  approximation  techniques  such  as  simulated 
annealing  heuristic  [90].  In  anticipation  of  finding  a more  global  minimum,  the 
simulated  annealing  heuristic  probilistically  accepts  configurations  which 
temporarily  deteriorate  the  quality  of  the  cost  function.  The  probability  of 
accepting  deteriorate  the  configurations  is  made  highest  at  the  early  iterations,  but 
is  reduced  slowly  as  the  iterations  progress,  thereby  accepting  far  fewer 
nonimproving  moves  near  the  minimum.  As  such,  a coarse  global  search  evolves 
into  a fine  local  search  near  the  global  minimum,  and  the  probabilistic  "jumps" 
provide  an  escape  from  local  minima. 

Here  is  the  analogy  between  the  annealing  simulation  and  the  optimal 
mapping  of  finite  elements  to  parallel  processors.  In  Eqn  (6.6)  where  the  cost  is 
a function  of  two  conflicting  factors  (computing  costs  and  communication  costs), 
in  the  annealing  problem  the  system  energy  is  a function  of  two  conflicting 
factors:  short  range  repulsive  forces  and  long  range  attractive  forces.  The  specific 
assignment  of  elements  to  processors  governs  the  total  cost  function  just  as  the 
atomic  positions  determine  the  body’s  energy  state.  In  the  hypothetical  situation 
in  which  the  attractive  forces  (communication  costs)  are  assumed  absent,  the 
system’s  energy  (cost  function)  may  be  minimized  by  the  new  dominant  repulsive 
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forces  (computing  costs)  requiring  that  the  atoms  (elements)  be  uniformly 
distributed  throughout  the  body  (among  processors).  Conversely,  if  the  repulsive 
forces  (computing  costs)  are  assumed  absent,  the  minimum  of  the  energy  (cost 
function)  will  be  achieved  by  placing  all  atoms  as  close  together  as  possible  (all 
elements  placed  in  one  processor).  Of  course  in  reality,  both  conflicting  factors 
are  present  and  their  effect  on  the  global  minimum  depends  on  their  relative 
strengths. 


6.3  Example  Problems 


As  an  example,  we  consider  a finite  element  representation  for  the 
vibrations  analysis  of  a laminated  composite  beam-plate  with  delamination.  The 
effects  of  in-plane  loading  on  the  vibration  characteristics  of  [02/902]s  and  [0/90]2s 
graphite-epoxy  composite  laminate  with  fixed-fixed  boundary  conditions  are 
examined.  For  simplicity,  the  delamination  is  assumed  to  be  located  symmetrically 
along  both  the  spanwise  direction  and  the  thickness  direction.  The  ply  mechanical 
properties  are  given  as  follows: 

E,  = 111.2  GPa,  E2  = 9.78GPa,  G,2  = 4.82  GPa,  v,2  = 0.32 
These  values  are  typical  of  an  advanced  fiber-reinforced  composite  materials. 
Similar  to  the  previous  chapters,  the  natural  frequency  as  well  as  the  in-plane 
loading  are  non-dimensionalized  in  the  following 
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where  I is  the  length  of  the  delaminated  plate,  m is  the  mass  per  unit  length  and 
D is  the  effective  bending  stiffness  of  the  laminate.  The  normalized  delamination 
length  is  defined  as  a/l  where  a is  the  delamination  length. 

Figure  6.2  shows  the  variation  of  the  fundamental  natural  frequency  with 
the  compressive  load  for  three  different  delamination  lengths.  It  can  be  clearly 
seen  that  the  fundamental  frequency  decreases  with  increasing  compressive  in- 
plane load.  A substantial  drop  of  the  frequency  becomes  more  substantial  at  the 
higher  compressive  load  although  the  effect  of  the  small  delamination  length  is 
negligible.  This  shows  the  effect  of  long  delamination  length  on  the  fundamental 
natural  frequency  is  important.  In  Figure  6.3,  the  effect  of  the  compressive  in- 
plane load  on  the  fundamental  frequency  with  two  different  layups  are  examined. 
It  shows  that  the  laminate  layups  also  have  some  effects  on  the  natural  frequency 


in  the  delamination  case. 
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Figure  6.2  Variation  of  Fundamental  Frequency  with  Compressive  In-plane  Load 

for  Different  Delamination  Lengths 


Figure  6.3  Variation  of  Fundamental  Frequency  with  Compressive  In-plane  Load 

for  Different  Different  Layups 


CHAPTER  7 

CONCLUSIONS  AND  FUTURE  WORK 


Refined  vibration  theories  of  sandwich  beam  and  plate  are  proposed  to 
investigate  dynamic  response  characteristics  of  sandwich  construction  with  newly 
developed  flexible  and  compressible  core  materials.  The  theory  is  employed  to 
study  the  vibration  damping  problems  of  sandwich  beam  and  plate.  The 
viscoelastic  face-core  interface  effects  on  structural  damping  of  sandwich  beam  are 
examined,  and  free  vibration  of  delaminated  sandwich  beam-plate  is  also  analyzed 
by  using  the  refined  theory.  Two  different  deformation 'modes  in  vibration  of 
delaminated  sandwich  beam-plate  are  discussed  and  the  physically  admissible  mode 
is  identified.  For  the  vibration  of  laminated  composites  with  delamination,  the 
layer-wise  finite  element  model  is  discussed  to  corporate  the  parallel  algorithm 
design  for  the  finite  element  analysis. 

7.1  Conclusions 

It  is  shown  that  the  proposed  theory  enhances  the  physical  insight  of  the 
vibration  damping  behavior  of  sandwich  beams  and  plates  consisting  of  more 
flexible  cores.  Parametric  study  demonstrates  that  the  shear  deformation  of  the 
viscoelastic  interface  layers  has  a significant  contribution  to  the  structural 
damping.  Effects  of  the  storage  shear  modulus  and  loss  factor  of  interface  layer 
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are  important  for  the  structural  damping  and  resonance  frequency.  Results  also 
demonstrate  that  the  there  may  exist  optimal  adhesive  bonding  stiffness  to  achieve 
the  maximum  sandwich  structural  damping.  This  theory  is  applicable  to  sandwich 
beams  and  plates  with  any  type  of  construction  and  flexible  core  materials. 

It  should  be  mentioned  that  as  damping  (loss  factor)  of  the  structure  is 
increased,  the  stiffness  of  the  structure,  in  general,  will  be  decreased.  Therefore, 
by  only  improving  damping  will  not  reduce  the  dynamic  responses  such  as 
displacements,  stresses  and  noise  level  of  the  structure.  It  is  a trade-off  situation 
between  damping  and  stiffness.  Generally,  a proper  value  of  interface  shear 
modulus  (adhesive  stiffness)  for  the  trade-off  between  enhancement  in  damping 
and  decrease  in  stiffness  can  be  obtained  in  most  situations. 

An  analytical  study  of  the  vibration  characteristics  of  sandwich  beams  with 
delamination  shows  that  the  effect  of  the  delamination  on  the  natural  frequencies 
depends  not  only  on  its  size  but  is  also  very  sensitive  to  its  spanwise  location  and 
vibration  mode.  Based  on  the  discussion  of  the  full  cycle  vibration  and  linear 
analysis,  it  can  be  seen  that  the  constrained  mode  model  is  a physically  feasible 
model  for  describing  the  vibration  behavior  of  delaminated  sandwich  beams 
whereas  the  free  mode  model  is  not.  The  layer-wise  finite  element  model  for  the 
laminated  composites  with  delamination  offers  the  possibility  of  dealing  with 
multiple  delaminations  and  complex  geometric  structural  problem,  in  addition  to 
accurate  analysis.  Its  expensive  computational  cost  brings  about  the  necessity  of 
implementing  parallel  algorithms  for  finite  element  analysis. 
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7.2  Future  Study  of  Related  Research 

A logic  extension  of  this  study  is  to  investigate  the  buckling  behavior  of 
delaminated  sandwich  composites  by  using  the  proposed  displacement  field.  Since 
the  energy  method  is  used  to  derive  the  governing  equations  of  motion  in  this 
study,  it  is  a consequence  to  formulate  the  energy  release  rate  equations  for  the 
fracture  analysis  of  delamination  growth  problems. 

Forced  vibration  problems  need  to  be  considered  for  the  vibration  damping 
of  delaminated  sandwich  beams  and  plates.  The  effects  of  delamination  on  the 
loss  factors  of  damped  sandwich  structures  are  required  to  address  in  the  future 
work,  and  the  correlated  experimental  study  must  be  performed.  It  is  noted  that 
there  exists  some  controversy  on  the  delamination  damage  effects  on  dynamic 
characteristics  of  laminated  composites,  thus  more  systematic  studies  in  this  area 
are  necessary. 

For  the  finite  element  analysis,  we  will  continue  to  study  the  multiple 
delamination  effects  of  laminated  composites.  The  study  can  extend  to  plate 
element,  and  more  study  on  parallel  algorithm  design  for  finite  element  analysis 
are  needed  to  meet  the  needs  of  growing  computational  demand. 


APPENDIX  A 

EQUATIONS  OF  MOTION  OF  A SANDWICH  BEAM 


As  we  mentioned  in  Chapter  2,  the  governing  equations  of  motion  and 
corresponding  boundary  conditions  of  a sandwich  beam  are  derived  via  the 
variational  principles  imposed  on  the  total  kinetic  energy  and  potential  energy  in 
terms  of  four  basic  functions,  w,,  W2,  Mq  The  governing  equations  of 

motion  can  be  written  in  the  following  matrix  form: 
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where  dMy  and  dKy  {i,  7 = 1,  2,  3,  4)  are  the  mass  matrix  differential  operators  and 


stiffness  matrix  differential  operators,  respectively.  This  equation  is  the  same  as 


Eqn  (2.23),  and  the  matrix  elemenets  are 
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The  consistent  boundary  conditions  are  also  derived  from  the  Hamilton 
principle  imposed  on  total  kinetic  and  potential  energy.  It  is  noted  that  the 
traction  boundary  conditions  and  mixed  boundary  conditions  given  in  Eqn  (2.24) 
through  (2.31)  in  Chapter  2 can  be  expressed  in  term  of  four  basic  functions, 
Wj(x,t),  W2(x,t),  Uo(x,t),  and  i|f(x,t).  All  boundary  conditions  at  two  ends  of  a 
sandwich  beam  are  given  as  follows, 
at  a:  = 0 and  x = I, 
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where  partial  differential  operators  dB]  and  dC|  {i  = 1,  2,  3,  4)  are 
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where  partial  differential  operators  dC')  (/  = 1,  2,  3,  4)  are 
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where  partial  differential  operators  and  dC]  (i  - 1,  2,  3,  4)  are 
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where  partial  differential  operators  aC*  (i  = 1,  2,  3,  4)  are 
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APPENDIX  B 

EQUATIONS  OF  MOTION  WITH  IMPERFECT  INTERFACE 


Since  the  imperfect  interface  conditions  were  considered  in  Chapter  3,  the 
Lagrange  multiplier  has  been  introduced  in  deriving  the  governing  equations  of 
motion  and  corresponding  boundary  conditions  of  a sandwich  beam  based  on  the 
variational  principles.  The  five  basic  functions,  W2(x,t),  u^{x,t),  U2(x,t)  and 

^(x,t),  are  used  in  the  governing  equations  of  motion  which  can  be  written  in  the 
following  matrix  form: 
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where  dMjj  and  dKy  {i,  7=1,  2. ..5)  are  the  mass  matrix  differential  operators  and 
stiffness  matrix  differential  operators,  respectively.  This  equation  is  the  same  as 
Eqn  (3.12),  and  the  elements  in  two  matrices  are 
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The  consistent  boundary  conditions  are  also  derived  from  the  Hamilton 
principle  imposed  on  total  kinetic  and  potential  energy.  Both  displacement  and 
traction  boundary  conditions  given  in  Eqns  (3.13)  through  (3-20)  can  be  expressed 
in  term  of  five  basic  functions,  M>2(x,t),  u^{x,t),  u^ix.t)  and  i|;(x,?),  as 

follows 

at  X = 0 and  x = I, 
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where  partial  differential  operators  dC]  (i  = 1,  2,. ..5)  are 
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where  partial  differential  operators  dB*  and  dC]  (i  = 1,  2,... 5)  are 
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